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' Abstract 

SPheno is a program that accurately calculates the supersymmetric particle spec- 
Qh. trum within a high scale theory, such as minimal supergravity, gauge mediated su- 

persymmetry breaking, anomaly mediated supersymmetry breaking, or string effective 
field theories. An interface exists for an easy implementation of other high scale mod- 
els. The program solves the renormalization group equations numerically to two-loop 
order with user-specified boundary conditions. The complete one-loop formulas for 
the masses are used which are supplemented by two-loop contributions in case of the 
neutral Higgs bosons and the [i parameter. The obtained masses and mixing matrices 
are used to calculate decay widths and branching ratios of supersymmetric particles 
as well as of Higgs bosons, b — * s'y, Ap and (g — 2) M . Moreover, the production cross 
sections of all supersymmetric particle as well as Higgs bosons at e + e~ colliders can be 
calculated including initial state radiation and longitudinal polarization of the incom- 
ing electrons/positrons. The program is structured such that it can easily be extend to 
include non-minimal models and/or complex parameters. 
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1 Introduction 

Supersymmetry (SUSY) [H |21 E] provides an attractive extension of the Standard Model 
(SM). It provides a qualitative understanding of various phenomena in particle physics: It 
stabilizes the gap between the Grand Unification scale / Planck scale and the electroweak 
scale jlj. It allows the unification of the three gauge couplings at a scale Mjj ~ 2 • 10 16 GeV 
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in a straight forward way jjp. The large top mass generates radiative electroweak symmetry 
breaking j^J. In addition it provides the lightest supersymmetric particle as a cold dark 
matter candidate jjj. Therefore, the search for supersymmetric particles is one of the main 
topics in the experimental program of present and future high energy colliders [HI El UD1 HI] • 

The Minimal Supersymmetric Standard Model (MSSM) consists of taking the Standard 
Model and adding the corresponding supersymmetric partners j^j. In addition a second 
Higgs doublet is needed to obtain an anomaly-free theory. The second Higgs doublet is 
also needed to give mass to u-type quarks and down-type quarks at the same time. The 
MSSM in its most general from contains more than 100 unknown parameters ^2] which are 
clearly to many for an exhaustive study. This number drastically reduces if one embeds the 
MSSM in a high scale theory, such as minimal supergravity theories jTHJ, gauge mediated 
supersymmetry breakingjH], or anomaly mediated supersymmetry breaking [15 . There is 
not yet a theoretical preferred scheme for supersymmetry breaking. For this reason it is 
important to know whether the precision of on-going and future experiments is sufficient: (i) 
To distinguish between the various schemes, (ii) To which extend it is possible to reconstruct 
the underlying theory. It has been demonstrated that the expected experimental accuracies at 
future e + e~ colliders complemented with data from the LHC allow for a successful reconstruct 
of such an underlying supersymmetric high scale theory ^H]- Connected with these questions 
is the question if the theoretical accuracy matches the experimental one. The present version 
of the program SPheno 1 is thought as a further step in getting accurate theoretical results to 
match finally the experimental precision. 

In the view of ongoing and future experiments it is highly desirable to have various and 
independent tools at hand performing the calculation of the supersymmetric spectrum, of 
decay widths, of branching ratios and of production cross sections. This allows for a cross 
check of the tools and by comparing the implemented methods and approximations one can 
also get a rough understanding of the theoretical error. SPheno is a program performing an 
accurate calculation of the supersymmetric spectrum, of the branching ratios of supersym- 
metric particles and the Higgs bosons and of the production cross sections of supersymmetric 
particles and the Higgs bosons in e + e~ annihilation including longitudinal beam polariza- 
tion. Moreover, the spectrum is used to calculate the branching of the rare decay b —>■ s'j, 
the supersymmetric contributions to the anomalous magnetic of the muon a M as well as 
supersymmetric contributions to the p parameter. 

For the calculation of the spectrum the programs ISA JET [T7], SOFTSUSY [TH] and 
SUSPECT [12] are widely used. A comparison of the results among these programs and with 
SPheno is given in [20] . The calculation of the branching ratios of supersymmetric particles as 
well as the production cross sections in e + e~ annihilation can be done with SPheno, ISAJET 
[T7j, SPYTHIA 21j, SUSYGEN and SDECAY j2SJ- A comparison of the results of these 
programs will be given in a future paper. 

SPheno has been written in Fortran90. The main focus has been on accuracy and on 
stable numerical results and less on speed. However, on a modern PC a typical running 
time is in the order of one second. The calculation is done using two-loop renormalization 
group equations (RGEs) [21j, complete one- loop correction to all SUSY and Higgs masses 
[23] supplemented by the 2-loop corrections to the neutral Higgs bosons [2H1 HI] and to 
the \x parameter |27j. The present version of SPheno does all calculations for real parameters 
neglecting the flavour structure in the fermion as well as in the sfermion sector. Decay widths 

1 SPheno stands for Supersymmetric Phenomenology 
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and cross sections are calculated using tree-level formulas. However, the couplings involved 
are running couplings and thus important numerical effects of higher order corrections are 
already taken into account. The program has been structured in such a way that the future 
inclusion of complex phases and mixing between the generations has already been considered 
in the design of the interfaces as well as in the definition of the various variables. Moreover, 
extensions of the MSSM, e.g. models with broken R-parity, can be implemented easily. 

The aim of this paper is to provide a manual of the program, version 2.0, to describe 
the approximations used and to display the results of a run. In Sect. [2] we will summarize 
the MSSM parameters and we give the tree-level formulas for the supersymmetric particles. 
Moreover, a short summary of the implemented high scale models is given. In Sect. El we 
list the implemented decay modes of supersymmetric particles and the Higgs bosons. We 
also discuss shortly the approximations used. In Sect. |U we present the implemented cross 
sections in e + e~ annihilation. In Sect. El we discuss give details on the implemented low 
energy constraints. In Sect. El we discuss the implemented algorithm in some detail. In 
Sect. |7|the main program is presented in detail providing the necessary information so that 
this program can be easily adapted to the user's requirement. In the appendices we dis- 
cuss the possible switches for influencing the program as well as a detailed discussion of 
possible input files. This can be done using SPheno specific files or by using the SUSY Les 
Houches Accord (SLHA) [22] • Moreover, we list the output of the program for a typical 
example. The source code as well as precompiled a version of the program can be ob- 
tained from the author via email: porod@physik.unizh.ch; or it can be downloaded from 
http : / /www-theorie . physik . unizh . ch/~porod/SPheno . html. 

2 MSSM parameters, particle spectrum, and Models 

In this section we fix our notation concerning the parameters and present the tree-level 
formulas for the masses as well as the mixing matrices. In the following we assume that the 
physical masses are ordered: m.j < rrij if % < j except for the sfermions as explained below. 
We also give a short overview over various high scale models which are implemented in the 
program. 

2.1 Ingredients for the Lagrangian 

The pure supersymmetric Lagrangian is specified by the Kahler potential giving the gauge 
interactions and by the Superpotential W giving the Yukawa interactions: 

W = e ab (YfL«HlE° + Y^Q^H^ + Y^Q b { H^ - »H?H») (1) 

where L, E, Q, D, and U denote the matter superfields. The SU{2)i representation indices 
are denoted by a, b = 1,2 and the generation indices by i,j = 1,2,3; e a b is the totally 
antisymmetric tensor with ex2 = 1. Note that the sign of \x is identical to the one in IS A JET 
[T7j and SOFTSUSY [TH] but opposite to the convention in [22]. Presently, real Yukawas 
Y L , Y D , Y u only are included. They and the gauge couplings are DR quantities. g\ is 
defined in the Grand Unification normalization g\ = Jh~/ 3g' where g' is the Standard Model 
hypercharge gauge coupling. 
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The next ingredient is the soft SUSY breaking Lagrangian, which is given by mass terms 
for the gauginos 

Lsoft,i = \ (M ± BB + M 2 W a W a + M 3 ~g a ~g a ) + h.c. , (2) 
mass terms for scalar matter fields and Higgs fields 

- M hAaQj ~ M hf*^ - M hA^ ( 3 ) 

and trilinear couplings of scalar matter fields and Higgs fields 

L sof t,3 = -e ab (A%L1E\E) + A%Q1H\D) + A%Q\H a 2 U] - B^H^ 

+h.c. (4) 



2.2 Masses and Mixing Matrices 

The masses of the various particles are induced by the soft SUSY breaking parameters and 
the vacuum expectation values V{ of the neutral Higgs fields Vi =< >. The ratio of the 
vacuum expectation values is denoted by tan/3 = v 2 /v\. The sum of the vacuum expectation 
values (vevs) squared is fixed by the gauge boson masses: 



2 1 2/ 2 , 2\ 

m w = 19 (v 1 +v 2 ), 



rrr z = \(gf i + <r)(v* + vZ 



(5) 



Neglecting the mixing between different generations, the Standard Model fermion masses are 
given by: 



1 



m. 



V2 



1 



V2 



1 
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(6) 



for u-quarks, d-quarks and leptons, respectively. 

The gluino mass is given by rrig = \M%\. The charginos are combination of the charged 
winos w ± = (w 1 =F iw 2 )/\/2 and the charged higgsinos h{ , h 2 . The Lagrangian contains the 
chargino mass term — (i/j~) t Xi/j + where ip~ = (—iw~,h^) T , ip + = (—iw + ,h 2 ) T and 



X = 



M ^v 2 



V2 



V 1 \1 



(7) 



The matrix is diagonalized by two unitary matrices U and V: 



(8) 



The neutral gauginos b, w 3 as well as the neutral higgsinos h®,h 2 form the neutralinos. In 
the basis ip° = (b , w 3 , h®, h 2 ) T one finds the mass term — (il)°) T Yip° with 



Y 



( M l 


i v 2 



o 

M 2 
9 2 v, 

-§1>2 



q' 

9 2 v, 


-n 



9 2 v 2 
-%v 2 
-fi 



\ 



(9) 



J 
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This matrix is diagonalized by an unitary matrix N: 

M D ^ = N*YN^ . (10) 

The CP-even electroweak eigenstates (H®, H®) are rotated by the angle a into the Higgs 
mass eigenstates (h°, H°) as follows: 



h \ _ I — sin a cos a \ I H± 
H° I ~ I cosa sina ) 1 H° 



(11) 



with m h o < m H o. The CP-odd and the charged Higgs masses are given by 

m 2 A o = B \x (tan (3 + cot f3) , m H+ — m^o + (12) 

at tree level. 

Neglecting generation mixing, the sneutrino masses are given by: 

Ml = Mf ij + im|cos2/3 (13) 
The other sfermion mass matrices are 2x2 matrices: 

IV1 !.i 1 (.. AL f..vL\*.. \ j/2 -2 _ __2 , —2 \ Lq: ) 



71 (viAfi - (fjY£)*v 2 ) M\.. - s 2 w c 2f3 m 2 z + m] 
M lu + (I " ¥w) <*P™z + << 75 (« 2 (^)* - »YV Vl 



1 7lM--Wf)^ 2 ) M^-I^ml + mlJ 1 J 

where c 2 {3 = cos 2/5 and sfy = sin 2 ^. These matrices are diagonalized by 2 x 2 matrices 
Rf ■ with 

mj- = RfMjR j f ~ (17) 

Sfermions are first ordered according to the generation and inside a generation according to 
their masses. For example, in the slepton sector the ordering is e±, e 2 , fj,±, jl 2 , f±, f 2 and 
similarly for squarks. 



2.3 High scale models 



In this section we summarize the key ingredients of the high scale models implemented in 
SPheno. We also present the formulas for the boundary conditions in the various models. In 
all cases the modulus is determined by requiring correct radiative symmetry breaking. At 
tree level the corresponding formula reads as: 



I |2 1 



tan 2(3 (m 2 H2 tan (3 - M 2 Hi cot (3) - m 



(18) 



Moreover, in all cases the high scale parameters are supplemented by the sign of \x and tan f3. 
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2.3.1 Minimal Supergravity 



The minimal supergravity (mSUGRA) scenario is characterized by a set of universal param- 
eters [T3J 12] at the GUT scale M GUT : the gaugino mass parameter M1/2, the scalar mass 
parameter M , and the trilinear coupling A : 

MiiMcuT) = M 1/2 (19) 
Mf(M GUT ) = M 2 (20) 
M m out) = A Y t (M GUT ) (21) 



2.3.2 Minimal Supergravity including right handed neutrinos 

In addition to the parameters of the mSUGRA model above the following parameters appear 
in this case: m UR , a common mass for all right handed neutrinos, and m Vi [i = 1,2,3), the 
light neutrino masses. In this case the MSSM RGEs are run up to the scale m UR where the 
neutrino Yukawa couplings are calculated using the formula Y u> i = y/m UR m u Jv2. In the range 
between m VR and Mqut the effect of neutrino Yukawa couplings is included in the RGEs of 
gauge and Yukawa couplings. At the GUT-scale the right sneutrino mass parameters as well 
as the trilinear coupling A Uji are given by: 

M&Mgut) = Ml (22) 
A v ,i(M GUT ) = A Y u>i {M GUT ) (23) 

The corresponding RGEs are used in the running from M G ut to m v R - At the scale m VR the 
neutrino Yukawa couplings Y„i, the trilinear couplings A v { and the soft masses M 2 Ri for the 
right sneutrinos are taken out of the RGEs and below the m VR the usual set of MSSM RGEs 
are used. 



2.3.3 Gauge Mediated Supersymmetry Breaking 

Gauge mediated supersymmetry breaking [T^ 13*0] (GMSB) is characterized by the mass 
Mm ~ (S) of the messenger fields and the mass scale A = (Fs) / (S) setting the size of the 
gaugino and scalar masses. The gaugino masses 

M 4 (M M ) = (N 5 + 3N w )g (A/M M ) a,(M M )A (24) 

are generated by loops of scalar and fermionic messenger component fields; Ni is the multi- 
plicity of messengers in the 5 + 5 and 10 + 10 vector-like multiplets, and 

1 + x 

g(x) = — — log(l + x) + (x — > —x) (25) 
x z 

is the messenger-scale threshold function 31 a which approaches unity for A Mm- Masses 
of the scalar fields in the visible sector are generated by 2- loop effects of gauge/gaugino and 
messenger fields: 

M|(M M ) = 2(N 5 + 3N w )f (A/Mm) E k l C i J a^(M M )A 2 (26) 

with ki = 1,1,3/5 for SU(3), SU(2), and U(l), respectively; the coefficients Cj are the 
quadratic Casimir invariants, being 4/3, 3/4, and Y 2 /4 for the fundamental representations j 



7 



in the groups i = SU(3), SU (2) and £7(1), with Y = 2(Q—I 3 ) denoting the usual hypercharge; 
also the threshold function [3T] 



/(«) = l + fflogd+.J-^^+ilJ.^)] 

+ (x -> -x) (27) 

approaches unity for A <C Mm- As evident from Eq. ("2""") scalar particles with identical 
Standard-Model charges squared have equal masses at the messenger scale Mm- In the 
minimal version of GMSB, the A parameters are generated at 3-loop level and they are 
practically zero at Mm- However, the program permits to set a value for A different from 
zero but universal for all sfermions. 



2.3.4 Anomaly Mediated Supersymmetry Breaking 

In anomaly mediated supersymmetry breaking (AMSB) the SUSY breaking is transmitted 
from the hidden sector to the visible sector via the super- Weyl anomaly ["""5] . The soft SUSY 
breaking parameters are given by: 

M a = -(3 a m 3/2 (28) 

9a 

A, = {3 Yi m 3/2 (29) 

M) = -7X/2 (30) 

where j3 a and /3yj are the beta functions of gauge and Yukawa couplings, respectively. The 
7j are the anomalous dimensions of the corresponding matter superfield and m 3 / 2 is the 
gravitino mass. Equation ("3~"")) leads to negative mass squared for the sleptons which is 
phenomenologically not acceptable. There are several possibilities to solve this problem [3*2*] 
and the simplest one is to add a common scalar mass M so that eq. (|3U|) reads as 

Mf = Ml^ 3 m\ l2 (31) 
This extension has been implemented in the program. 



2.3.5 String Effective Field Theories 

Four-dimensional strings naturally give rise to a minimal set of fields for inducing supersym- 
metry breaking; they play the role of the fields in the hidden sectors: the dilaton S and the 
moduli T m chiral superfields which are generically present in large classes of 4-dimensional 
heterotic string theories. The vacuum expectation values of S and T m , generated by genuinely 
non-perturbative effects, determine the soft supersymmetry breaking parameters [3*3*1 134| . 

In the following we assume that all moduli fields get the same vacuum expectation values 
and that they couple in the same way to matter fields. Therefore, we omit the index m 
and take only one moduli field T. The properties of the supersymmetric theories are quite 
different for dilaton and moduli dominated scenarios as discussed in The mass scale 

of the supersymmetry parameters is set by the gravitino mass 7713/2. 

In the program we implemented the complete 1-loop formulas given in [31] . Three classes 
of models are implemented in the program: two versions of Oil compactification defined by 
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the sets A and B of boundary conditions in [34] as well as an 01 compactification scheme. 
For the implementation of the 01 compactification scheme we have used formulas Eqs. (3.21) 
- (3.23) of 



Mi 



-9i m 3/2 



VSsm9 + 



bi + s V3 sin 0g 2 s Ci-J^C? 



M? 

3 



+ 2tcos9G 2 (t) [S GS + b i -2 y £C{(l 



/167T 2 



(32) 



2s 



3/2 1 (l + n j c °s 2 9) + sin 9 

+ 7i + 2* cos 9 G 2 (t) lj m ( n J + n k + n m + 3)\ 

km ) 



A; 



jkm 



m 3/2 



V3sm9 - 2tcos6(nj + n k + n m + 3)G 2 (t) + jj + 7fc + 



(33) 
(34) 



s = (S) is the vacuum expectation values of the dilaton field, t = (T)/m 3 / 2 is the vacuum 
expectation value of the moduli field(s), and G%(t) = 2((t) + l/2t is the non-holomorphic 
Eisenstein function with ( denoting the Riemann zeta function. 5gs is the parameter of the 
Green-Schwarz counterterm. jj are the anomalous dimensions of the matter fields, the 7] 
and 7^ m are their gauge and Yukawa parts, respectively. Cj, C{ are the quadratic Casimir 
operators for the gauge group Gi, respectively, in the adjoint representation and in the matter 
representation. The indices i,j, k denote Hi, H2, E, L, D, U and Q. The A-parameters are 
finally given by: 



A e>n (GUT) = Y e , nn (GUT)Az nLnHi 
Ad, n (GUT) = Y d , nn (GUT)A 6n Q nHi 
A u , n (GUT) = Y U:nn (GUT)A &n Q nH2 



(35) 
(36) 
(37) 



where n denotes the generation. 

In case of the Oil compactification scheme the gaugino masses are given by Eqs. (3.11) 
of HI]: 



Mi 



-9i m 3/2 



jv^sinfl 



2k 



1/2 



+ 



16tt 5 



2t cos 9G 2 {S G s + bi) + bi 
V3g 2 s sin 9 



2k 



1/2 



(3* 



For the sfermion parameters we have implemented two sets of boundary conditions: set (A) 
is specified by formulas Eqs. (3.15) and (3.19) of 



M 



m 



3/2 



sin 2 9 + 7i + 



\^3sm9 



k„-r 



1/2 



A 



"i3/2S7i + 7i + 7fc- 



sin# 



k 



1/2 



(39) 
(40) 
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Set (B) is specified by formulas Eqs. (3.16) and (3.20) of 



M 



ni 



3/2 



f \/3sinfl 

+ sin 2 
9 sin 2 9 



1 + 2t cos 6G 2 



z jk 



1 + 7j + In [2*|t?(*)| 4 ] ( E 7f + 2 E 7f ) - E 7? M<7a) 

a j'fc a 



+ In [(t + t)\ V (t) | 4 ] ( E 7f (t) + 5 E 

4 J # 



.4 



"^3/2 1 (7. + 7j + 7fc) 
v^sin^ 



k 



1/2 



1 + 2* cos#G 2 

2 



^ + E|(7r + 7j + 7 fc a )(l-ln(^)) 



-In 



(t + t)|r7(t)r 



(E^(7? + 7* + 7 fc a ) 

a 

EM7i m + 7] m + 7l m )) 



Zm 



(41) 



(42) 



In all three cases we have assumed that terms proportional to the log(/2j) can be neglected 
(p,i denote the Pauli Villar masses). 



2.3.6 General High Scale Model 

It is clear from the examples above that up to now there is no unique mechanism for super- 
symmetry breaking. Therefore, we have implemented the possibility to specify rather freely 
a high scale model. This model is specified by: a set of three in principal non-universal 
gaugino mass parameters M 1 / 2 [U(1)}, M 1 / 2 [SU(2)}, M 1 / 2 [SU(3)}; a scalar mass for each type 
of sfermion, resulting in fifteen parameters: M% .., M~ .., M% .., M% .., M% . .; two Higgs mass 

parameters Mjj and Mjy ; nine different A parameters A),e,iij A ^,u and A 0)U> a. Here ii 
denotes that only the diagonal entries can be set, because in the current version the effects 
of generation mixing is not taken into account. A model of this kind has been used in 
for the study of low energy observables and the supersymmetric spectrum. It also can be 
used, for example, to set the boundary conditions for the gaugino mediated supersymmetry 
breaking [3^]. This general model will be denoted by SUGRA. 



2.3.7 General MSSM at low energies 

Starting with version 2.2.0 there exists also the possibility to give the parameters at the 
low scale M EWSB together with scale. In this case the parameters are taken as running 
parameters at the scale M EWS b and the masses and mixing angles are calculated using these 
parameters which in turn serve as input for the calculation of decay widths and cross sections. 
The input parameter are: the electroweak gaugino mass parameters Mi and M 2 , the gluino 
mass rrif the parameters describing the Higgs sector /i, tan/3, the mass of the pseudoscalar 
Higgs mass m A ; the sfermion mass parameters M Ejii , M~ L ii , M b ii , Mfj u , Mq u , A^, A dtii , 
and Aiji. This general model will be denoted by MSSM. 
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3 Decays of supersymmetric particles and Higgs bosons 



The programs calculates the most important two- and three-body decays of supersymmetric 
particles at tree level. In case of three-body decays the formulas are implemented such, that 
the effects of decay widths in the propagators are taken into account [Slj. Therefore, it is 
possible to perform the calculation even in case that some of the intermediate particles are 
on-shell. This is useful in the case that the two-body decays have small phase space, because 
then the calculation of the three-body decays gives a more accurate result, e.g. T(xt —> 
X° 1 W + )xBR(W + — > vl + ) can be quite different from T(xf —> Xivl + ) if the decay xt —> xV^ + 
has only small phase space. 

The following sfermion decays are calculated: 

ft - fxl fxt (43) 
/, - ^ Z°, /; W ± (44) 

U -> ] 3 (h\R\A»), f' J W ± (45) 

In case of the lighter stop, it is possible that all two-body decays modes are kinematically 
forbidden at tree-level. In this case the following decay modes are important |3*H1 I39| I4*U]: 

h -> cx?, 2 (46) 
h - W + bxl H + bx\ (47) 
t x -> buff, bl + D (48) 

where I — e, fi, r. The corresponding widths are calculated within SPheno using the formulas 
given in [39| . In case of GMSB models scenarios exist where the charged sleptons are next 
to lightest supersymmetric particles (NLSP) and the gravitino G is the LSP. In this case the 
sleptons decay according to: 

k -> IG (49) 

Here we use the formulas given in [30J. 

It is well known that the partial widths of sfermions can receive considerable radiative 
corrections [H] . However, the branching ratios are not that strongly affected [12] • Therefore, 
for the moment being tree-level formulas are implemented. Some important numerical effects 
of higher order corrections are nevertheless implemented by using 1-loop corrected masses and 
running couplings in the formulas. The complete implementation of higher-order corrections 
is left for future versions of the program. 

In case of charginos and neutralinos the following decay modes are calculated: 





- Z x p W ± xl 




(50) 


"K.% 


-> {h\H\A»)xl 




(51) 




~* ffji ffj 




(52) 


xt - 


-> z°xt, W + x] 




(53) 


xt - 


-> (h°, H°, A ) xt-. 


H+x] 


(54) 


xt - 


- ff'i 




(55) 
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In case that all two body decay modes are kinematically forbidden the following three-body 
decays are calculated: 



X* - ffxl ffxl (56) 



x° -> qqg (57) 

Xt - ffxt, ff'x" (58) 
Xk ~> 93 5 (59) 

In the calculation we have included all contributions from gauge bosons, sfermions and Higgs 
bosons [37J ESI- The Higgs contributions can be important in certain regions of parameter 
space jH]. In addition the loop induced decay s 

X- - X°7 (60) 

are calculated jlH] taking into account the left-right mixing of sfermions. Similarly to case of 
the sleptons there exist parameter regions in GMSB models where the lightest neutralino is 
the NLSP and it decays according to 

X? - 7 G (61) 
Xi ^ Z° G (62) 
X? - h° G (63) 

Here we use the formulas given in [30j. 

In case of gluinos the following two-body decays are calculated: 

9 -»• qq% (64) 

with q = u,d, c, s, t, b. Again, in case that these decays are kinematically suppressed, the 
three-body decay modes are calculated: 

9 X°9<? (65) 

9 xfi'q (66) 

g -> bW't u bW + tl (67) 

Here we have implemented the formulas given in [4*6] . In addition the decays 

^ ^ X-7 (68) 

are calculated |431 117] taking into account the left-right mixing of sfermions. 
In case of Higgs bosons the following decays are calculated: 

(69) 
(70) 
(71) 
(72) 
(73) 
(74) 



4> - 


- // 


<t> - 


~* fi f j 






<t> - 


-> xlxl 


<t> - 


~* "Kr "Ks 


H° - 


-> z°z°, w + w 


H° - 


-> h°h° 
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A -> h°Z° (75) 

- //' (76) 

# + - /i/; (77) 

# + - X^x; (78) 

# + -> (79) 

with = ft, , if , A and / = z^, e, /i, r, w, d, c, s, t, 6. It is well known, that the widths as well 
as the branching ratios of the Higgs bosons can receive large one-loop corrections [JBl EE ED] ■ 
In the present version only the gluonic QCD corrections for the decays into quarks [IB] have 
been implemented. Therefore, the numbers provided by SPheno have to be taken with care 
and for refined analysis other programs, such as HDECAY 51J should be used. 

4 Production of supersymmetric particles and Higgs 
bosons 

The program calculates the following cross sections: 

e + e- -> fjj (f = l,v,q) (80) 

e + e- - xlxl (81) 

e + e- - x+Xs (82) 

e+e" -> h°Z, H°Z (83) 

e + e- -> (84) 

e+e" (85) 

We haven taken the formulas of for sfermion production, jSHl EH] for production of 
charginos and neutralinos and (SD] for Higgs boson production. Initial state radiation has 
been included using the formula given in In case of squarks in addition QCD corrections 
due to gluon exchange are included |HH EHj. Care has to be taken in case one calculates 
the cross sections near threshold because then higher order corrections are important to get 
reliable results and, thus, the numbers obtained in the program have to be taken with 
care near the threshold. All cross sections are implemented such, that one can specify the 
degree of longitudinal polarization P e - of the incoming electron beam as well as the degree 
of longitudinal polarization P e + of the incoming positron beam. Here P e - is within the range 
[—1, 1], where {—1, 0, 1} denote 100% left-handed electrons, completely unpolarized electrons 
and 100% right-handed electrons, respectively. The same notation is used in case of positrons. 
For example, P e - = —0.8 (P e + = —0.8) means that 80% of the electrons (positrons) are left- 
polarized whereas the remaining 20% are unpolarized. 



5 Low Energy Constraints 

The supersymmetric parameters are constrained by direct searches at colliders and by loop- 
effects which supersymmetric particles induce observables of low energy experiments. Pro- 
vided one neglects mixing between different sfermion generations the following quantities 
constrain several parameters of the MSSM: the rare decay b — > s 7, the anomalous magnetic 
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moment of the muon a M and the supersymmetric contributions to the p parameter. These 
constraints are implemented in the program using the formulas given in jJ7| EH] for b — > s 7 
supplemented by the QCD corrections as given in jSH], |)U] for a M and [HI] for the sfermion 
contributions to the p parameter. In call cases we use the running couplings at rriz for the 
calculation of the observables. The use of running couplings together with the correct im- 
plementation of supersymmetric threshold corrections for the couplings results in taking into 
account the most important higher oder corrections as has been pointed out e.g. in EH] 
for the case of b — > sj. The implementation of the supersymmetric threshold corrections to 
the couplings will be discussed in the next section. 

6 Details of the Calculation 

In this section we describe the algorithm used. It is schematically displayed in Fig. ^ The 
following standard model parameters are used as input: fermion masses, the Z-boson pole 
mass, the fine structure constant a, the Fermi constant Gf and the strong coupling constant 
a s (mz)- It is assumed that a s {mz) is given in the MS scheme. We describe first the 
implementation of the high scale models and comment then on the case of the implemented 
MSSM model. 

6.1 First rough calculation of SUSY and Higgs boson masses 

In a first step, we calculate gauge and Yukawa couplings at m z scale using tree-level formulas. 
These are used as input for the one-loop RGEs to get the gauge and Yukawa couplings at 
the high scale where also the boundary conditions for the high scale model under study are 
imposed. Afterward one-loop RGEs are used to get a first set of parameters at the electroweak 
scale. These parameters are used to get a first set of supersymmetric particle masses and 
Higgs masses using tree-level formulas except for the neutral CP-even Higgs bosons where 
one-loop effects due to (s)quarks of the third generation are taken into account. These masses 
are the starting point for the iterative loop which calculates the spectrum within the required 
precision as described below. 

6.2 Main loop for the calculation of SUSY and Higgs boson masses 

In the next step the gauge couplings and sin 2 6 W are calculated at m| in the DR scheme 
using the formulas given Appendix C of [25]. The only difference is that we use an updated 
number for Aai ep + Aahad = 0.06907 [62] . The mass of the W-boson is calculated using the 
formula [HH] 



ml 



w 




(86) 



1 



(87) 



P 



1 - Ap 




(89) 



8. 



.8 
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Calculate gi(M z ), h tjbjT (M z ) at tree level. 



I 



Run to M x at 1-loop. Apply soft SUSY breaking boundary condition 



I 



Run to M EWS b- Calculate and sparticle pole masses at tree level 



SM and SUSY radiative corrections to gi(M z ), h tt b,r(M z ) 



I 



Run to M X - Apply soft SUSY breaking boundary condition 



Run to M EWSB . Calculate and sparticle pole masses at loop level 



I 



Check if required precision is achieved 



no 



yes 



Calculate decay widths and branching ratios 



Calculate cross sections 



I 



Calculate low energy observables 



Figure 1: Algorithm used to calculate the SUSY spectrum, decay widths, production cross 
sections and low energy observables. Each step (represented by a box) is explained in the text. 
The initial step is the uppermost one. M EWSB is the scale at which the EWSB conditions are 
imposed, and M x is the scale at which the high energy SUSY breaking boundary conditions 
are imposed. 
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where Hyy is the transverse part of the vector boson self-energy and 5vb contains non- 
universal corrections to the /i decay. In the computation we have included the leading SM 
two- loop corrections the formulas of and adapting the Higgs contribution as in [25] . The 
formulas of the SUSY contribution to 5vb are given in [U71 [HE] • All masses appearing in the 
loops are running except for the top-quark, because the 2-loop part is given for an on-shell 
definition of the top mass. Note here that also the gauge boson masses in the loops are 
running masses and thus an iteration has to be performed in practice. 

For the calculation of the Yukawa couplings we use the complete formulas for the fermion 
masses and the vacuum expectation values given in Appendix D of [25]- In addition we have 
implemented the following improvements. The five light quarks and a s are evolved to mz 
using 4- loop RGEs as given in {6H| including threshold corrections at the low scale [ZOj- This 
evolution is done in the MS scheme. At mz the shift to the DR scheme is performed using 
the formulas given in [7TJ[72], e.g.: 

m vm {m z ) = m hM - s (m z ){l-^- |^ + ^ - (90) 

where a s is given in the DR scheme which is the reason for the different factor in front of a 2 s 
compared to [12] . We use the complete formulas given in Appendix D of [25. to calculate the 
SUSY contribution denoted by Arrib t susY(mz)- For the calculation we use running gauge 
and Yukawa couplings at mz- The SUSY contributions Am^suSY^z) are then resumed 
using [73] 

m bm {m z ) = - AmbSUSY(mz) (91) 

m b.DR (mZ ' ) 

In the case of the top-quark a modified procedure has to be used due to its large mass. 
In this case we start from the pole mass m< and use the formulas given in [25] to obtain 
mf R (mz). The difference compared to [25] is in the 2-loop parts where we have taken the 
exact 2-loop formula given in [7Tj : 

= Y^WWsr 2011 16A2 _ 803) + 246 \ 
* V^/V 918 9 3 3 J K ' 

L = log (4) ( 9 3) 

In the case of leptons we first calculate the running masses mi(m{) at 1-loop level in the 
MS scheme. Afterward we use two-loop RGEs as given in [7j] to evolve them to mz- Here 
we perform the shift to the DR scheme using the formula: 

DR,SM I \ ~ms,sm, \ 3(g /2 — 
m { {m z ) = m l (m z ) II I (94) 

The analog of Eq. (|91|) is used to obtain the running mass including the SUSY threshold 
corrections. 

The obtained gauge and Yukawa couplings are evolved to the high scale using two-loop 
RGEs |24j . The high scale can either be fixed or can be calculated from the requirement 
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9\ = 92 at the high scale. The various cases are discussed in Sect. [7|and Appendix |XJ At 
the high scale the boundary conditions for the soft SUSY breaking parameters are set. The 
implemented models are summarized in Sect. 12.31 see also Appendix IB. 31 The complete 
set of two-loop RGEs 24 is used to evolve the parameters down to the electroweak scale 
Mew sb — ■\/ m ti m h (there exists also the possibility to fix M EWS b to a constant value as 
described in Appendix IA1 . 

The parameters are used as input to calculate the sparticle pole masses at one-loop order 
and in case of the neutral Higgs at two-loop order. Here we use the complete formulae given in 
the appendices of for the one-loop contributions and for the 2- loop corrections 0(a s (a t + 
ab) + (at+ab) 2 +a T (yb+Oifau) for the neutral Higgs boson the formulas given in [SEIEZIEH]- For 
the 0(a s otb) contributions we use the complete expressions which can be obtained from the 
0(a s a t ) contributions by appropriate replacements. In case of sfermions we have included in 
all cases left-right mixing. All gauge and Yukawa couplings are understood as DR quantities 
at M EW sb- Also tan/? and the vacuum expectation values are evaluated at M EWS b to get a 
consistent set of input parameters. Note, that we express in all couplings the fermion masses 
and gauge boson masses by their corresponding expressions due to gauge couplings, Yukawa 
couplings and vacuum expectations values in the formulas of |25j . In all cases running masses 
are used as input for the loop integrals. In addition we have implemented the 0(a s (a t + 
ctb) + (a t + a>b) 2 + a T ab + afau) corrections for the calculation of [23I2H1- The numerical 
evolution of the one-loop integrals is based on the FF package and the Loop Tools package 

m. 

The masses obtained are used as input for the next iteration which starts again by calcu- 
lating the SUSY contributions to gauge and Yukawa couplings at ra z . In the case that during 
this iterative process an unphysical situation occurs, e.g. a pole mass squared being negative, 
the program terminates and it provides information on the exact reason for termination. The 
iteration is continued until all relative differences between the sparticle masses are smaller 
then the user imposed quantity 5: 

\rrii — , . 

5 > l — i-^i 95 

rrii 

for all sparticle masses; i denotes the i-th iteration. In most cases this achieved after three to 
four iterations. In the case that more than the maximal allowed number of iterations (user 
specified) are necessary, the program leaves the iteration giving a warning message. 

In the case of the MSSM model the running between M ew $b and M x is obviously omitted. 
In this model one has, however, to preform an iteration to get the Yukawa and gauge couplings 
correct because tan/3 is defined at M E wsb instead of mz- 



6.3 Calculation of the other observables 

The masses and mixing angles are then used to calculate the branching ratios and decay 
widths. Here two- and three body decays of supersymmetric particles are calculate. Note 
that we use the couplings as input which are renormalized at the scale M EWSB . The user 
has the possibility to force the program to calculate three body decays even if one or more of 
the intermediate particles are on-shell. This possibility is useful in the case where the 2-body 
decay has only small phase because then the calculation of the three body decay width(s) 
give a more reliable result. 
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Afterward the production cross sections at an e + e~ collider of all supersymmetric particles 
as well as all Higgs bosons are calculated. Here the user has the possibility to specify the 
center of mass energy as well as the degree of longitudinal polarization of the incoming beams. 
Moreover, the user can specify if initial state radiation shall be included in the calculation or 
not. 

Finally, the low-energy constraints described in Section are calculated: BR(b — > sj), 
SUSY contributions to a M and the sfermion contributions to Ap. For the the calculation 
of these quantities we evolve the gauge and Yukawa couplings from the scale M EWSB down 
to mz- The couplings at rriz are then used as input for the calculation of the low energy 
observables. For example in calculation of BR(b — > sj) the most important contributions to 
the Cf coefficients are implemented as 



K ts KtbXtw /2 
4< 



KtsKtb 
4m 2 H+ 



3 

Y t 2 cos 2 p /2 
4 V3 



Fi(xtw) + F 2 (x t w) 

Fi(xtH+) + F 2 (x tH +) 



(96) 



Y b Y t cos P sin (5m t f 2 . . 

o F z( x tH+) + F A (x tH i 

m b \3 



(97) 



- K ts K tb 



y 



-C L! ijC R>ij ( -i^O^tt.) + ^(^ 



YbRt,aUj2 

-gRt^Vji + Y t R~ m V j2 



(98) 

(99) 
(100) 



Here Yj are the Yukawa couplings, U and V are the chargino matrices, K is the CKM 
matrix, Ri is the stop mixing matrix and x ab = m 2 a /m 2 . The loop functions Fj are given in 
|57j . A similar replacement is done in the contributions to the Cg coefficient. Moreover, in 
the program also the contributions from the first two generation of (s)fermions to are 
included for completeness. We then use j^H] to obtain 



BR{b -> S7) 



1.258 + 0.382r7 2 + 0.015r 2 



1.395r 7 + 0.161r 8 + 0.083r 7 r 8 



(101) 



where r 7 = Cj/Cj(W + ) and r 8 = Cs/Cs(W + ). In this way important higher order corrections 
are taken into account, in particular the large tan/? effects in case of b — > 37 



7 A sample example 

In this section we discuss the executable statements of the main program given in the file 
SPheno . f 90. In the first statements the required modules are loaded and the various variables 
are defined. Afterward the error system is initialized and the input data are read by calling: 

Call ReadingData(kont) 

The source code can be found in the file SPheno. f 90. This routine checks first if the file 
LesHouches . in exists where the input data are provided according to SLHA (for a detailed 
description see Appendix ■ If this file is absent, the routine checks if the files Control . in 
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Table 1: Variables for parameters and couplings. The parameters are defined in the module 
MSSM_Data and they are explained in Sectional dp means double precision, 
parameter/coupling type & Fortran name 



complex (dp) : : phase_mu 

tan/3 real (dp) : : tanb 

M h M 2 ,M 3 complex(dp) :: Mi(3) 

Mf, Ml complex(dp), dimension(3,3) :: M2_E, M2_L 

M|,, Mq, My complex(dp), dimension(3,3) :: M2_D, M2_Q, M2_U 

Ai,Ad,A u complex (dp) , dimension(3,3) :: A_l, A_d, A_u 

fi complex (dp) : : mu 

Bfi complex(dp) : : B 

M 2 H real (dp) : : M2_H(2) 

g', g real (dp) : : gp, g 

Yi,Ya,Y u complex(dp), dimension(3,3) :: Y_l, Y_d, Y_u 

v 1 ,v 2 real(dp) :: vevSM(2) 

g',g,g s real (dp) :: gauge (3) 



and/or StandardModel . in are present. The first one can be used to set various flags (see 
Appendix IB.1|) whereas the second one is used to set the SM input (see Appendix IB.4|) . 
Standard values as described in the Appendix are used if any of these files is absent. Afterward 
the model specific information is read from the file HighScale . in which is described in 
Appendix IB. 31 (for a short description of the implemented models see Sect. 12. 3|) . In all 
cases the file Messages. out is opened at channel 10 where all warnings and/or debugging 
informations are stored. 

Before calling the subroutine CalculateSpectrum, which performs the calculation of the 
spectrum, the user has the possibility to fix the high scale and/or the scale where the pa- 
rameters and the loop corrected masses are calculated. For this purpose one or both of the 
following lines must be uncommented in the program: 

! Call SetGUTScale(2.el6_dp) ! please put the GUT scale 

! Call SetRGEScale(l.e3_dp**2) ! please put the scale M.EWSB squared 

The default is that these scales are calculated by the program. The high scale is computed 
from the requirement g\ = g 2 (except in GMSB where the high scale is an input). The scale 
Mewsb is given by M E wsb = ^m {l m i2 . 

The accurate calculation of the SUSY parameters and the spectrum is done by the fol- 
lowing call: 

delta = l.e-3_dp 
WriteOut = .False. 
n_run = 20 

If (kont.Eq.O) Call CalculateSpectrum(n_run, delta, WriteOut, kont, tanb, vevSM & 
& , mC, U, V, mN, N, mSO, mS02, RSO, mPO, mP02, RPO, mSpm, mSpm2, RSpm & 
& , mSdown, mSdown2, RSdown, mSup, mSup2, RSup, mSlepton, mSlepton2 & 
& , RSlepton, mSneut, mSneut2, RSneut, mGlu, PhaseGlu, gauge, Y_l, Y_d & 
& , Y_u, Mi, A_l, A_d, A_u, M2_E, M2_L, M2_D, M2_Q, M2_U, M2_H, mu, B & 
& , m.GUT) 
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Table 2: Variables for masses and mixing matrices as given by the routine 
CalculateSpectrum. They are defined in the module MSSM_Data and their connection to 
the parameters at tree-level is explained explained in Section |21 dp means double precision, 
masses / mixing matrix type & Fortran name 



m- g 

u,v 

N 

m h o, 
m G o, 

R/3 

m G + 
R'p 

VTiy 

Ru 

m l 
R-i 

ma 



rriHO 



rriAo 



Ru 
Ri 



real (dp) :: mglu 
complex(dp) :: PhaseGlu 
real(dp) :: mC(2) 
complex(dp) :: U(2,2), V(2,2) 
real(dp) :: N(4) 
complex(dp) :: N(4,4) 
real(dp) :: mS0(2) 
real(dp) :: RS0(2,2) 
real(dp) :: mP0(2) 
real(dp) :: RP0(2,2) 
real(dp) :: mSpm(2) 
complex(dp) :: RSpm(2,2) 
real (dp) :: mSneut(3) 
complex(dp) :: Rsneut(3,3) 
real(dp) :: mSlepton(6) 
complex(dp) :: Rslepton(6,6) 
real(dp) :: mSup(6) 
complex(dp) :: Rsup(6,6) 
real(dp) :: mSdown(6) 
complex(dp) :: Rsdown(6,6) 



The meaning of the various variables and their type is given in Tables ^ and 121 Variable 
names ending with "2" indicate masses squared. The variables for the mixing matrices are 
already structured for a latter extension to include the effects of generation mixing and/or 
complex phases: the sfermion mixing matrices are 6x6 (except for sneutrinos which is a 
3x3 matrix) . In the present release most of the entries are zero except for the diagonal 2x2 
blocks which contain the left-right mixing for every species of sfermions. For example, the 11, 
12, 21 and 22 entries in Rslepton specify the left-right mixing of selectrons, and similarly 
the 33, 34, 43 and 44 (55, 56, 65 and 66) entries specify the left-right mixing of smuons 
(staus). The squark mixing matrices have the same generation structure. This structure has 
been chosen to facilitate a later extension which includes flavour violating entries. Beside the 
variables given in Tables ^ and El the following variables appear: 

• delta : specifies the required relative precision on the masses. If the maximal relative 
differences between the physical masses obtained between two runs is smaller than 
delta, the routine Sugra leaves the iteration loop. 

• m_GUT : the value of the scale where the high energy boundary conditions are imposed. 

• kont : A variable which is provided everything is o.k. Otherwise either a numerical 
problem has occurred and/or the parameters belong to an unphysical region, e.g. a 
minimum of the potential where charge and/or colour breaking minima occur. In such 
a case the information is written to the file Messages . out. 
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• WriteOut : If it is set .True, then intermediate debugging information is written to 
the screen and the file Messages . out. 

• n_run specifies the maximal number of iterations of the main loop. A warning will be 
given in the case that the required precision delta has not been reached within n_run 
iterations. 

Note that the parameters are running parameters at the scale Mewsb- The complete spec- 
trum is calculated at 1-loop level using the formulas given in [22]. The exceptions are the 
masses of the neutral Higgs bosons (scalar and pseudo-scalar) as well as \i the two loop 
corrections due to a s and all third generation Yukawa couplings are included [26| I27j. 

In the next part the branching ratios, the partial decay widths and the total decay widths 
are calculated provided that L_BR= . TRUE . and kont.eq.O: 

If ((L_BR) .and. (kont.eq.O)) then 
epsl = l.e-5_dp 
deltaM = l.e-3_dp 
CalcTBD = .False. 
ratioWoM = l.e-4_dp 

If (HighScaleModel . Eq. "GMSB") Is.GMSB = .True. 

Call CalculateBR (gauge, mGlu, PhaseGlu, mC, U, V, mN, N & 



& 


mSneut, RSneut, mSlepton, RSlepton, mSup, RSup 


& 


& 


mSdown, RSdown, uL_L, uL_R, uD_L, uD_R, uU_L, uU_R 


& 


& 


mSO, RSO, mPO, RPO, mSpm, RSpm, epsl, deltaM 


& 


& 


CalcTBD, kont, ratioWoM, Y_d, A_d, Y_l, A_l, Y_u 


& 


& 


A_u, mu, vevSM, Fgmsb, m32 


& 


& 


gP.Sl, gT_Sl, BR_S1, gP_Sn, gT_Sn, BR_Sn, gP_Sd 


& 


& 


gT_Sd, BR_Sd, gP_Su, gT_Su, BR_Su, gP_C, gT_C 


& 


& 


BR_C, gP_N, gT_N, BR_N, gP_G, gT_G, BR_G, gP_P0 


& 


& 


gT.PO, BR_P0, gP_S0, gT_S0, BR_S0, gP.Spm, gT.Spm 


& 


& 


BR.Spm) 





end if 



Variables starting with gP_, gT_ and BR_ indicate partial widths, total widths and branching 
ratios, respectively; they are Real (dp) vectors. The first index is the index of the decaying 
particle whereas the second one gives the mode. The correspondence between the second 
index and the modes is summarized for sfermions (variables ending SI, Sn, Sd and Su for 
sleptons, sneutrino, d-squarks and u-squarks, respectively) in TableHH for charginos in TableEJ 
for neutralinos in Table for gluinos in Tableland for the Higgs bosons in Tables [7| and |H1 
Here the following variables are new: 

• epsl : gives the accuracy to which the phase space integrals in three body decays are 
calculated. 

• deltaM : this variable affects the calculation of the phase space integrals in three body 
decays. In case that mi/ (m — J2i mi) < deltaM than m, is set to zero in the calculation 
of the phase space integrals, m denotes here mass of the decaying particle and 
(i=l,2,3) are the masses of the decay products. 
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Table 3: Correspondence between the indices for sfermion partial widths (branching ratios) 
and the modes. t\ is listed extra because it can have three-body decay modes. 



mode 


Z 


V 


d 


u 




fi -> flft 

J i J A. k 


1-4 


1-4 


1-4 


1-4 


1-4 


fi -> f'Xj 


5-6 


5-6 


5-6 


5-6 


5-6 


fi -> fa 






7 


7 


7 


f. _> w± f ' 

f . ZT± f / 

h H h 


7 


7-8 


8-Q 




8-Q 


c 
o 


q in 


1 n 1 1 

1U- 1 1 


1 n 1 1 

1U- 1 1 


10 11 
1U-1 1 


J 2 — >■ ^ Jl 


Q 










r /I r 


10 




13 


13 






11 


- 


14 


14 


- 


/ 2 -> #°A 


12 




15 


15 


- 


J 7 /*~Y 

<i — > / G 


13 










*i C X°2 










55-56 


fi -> W + 6x? 










57 


fi — > 6e + z> e 










58 












59 


ii — > bT + v T 










60 


fi -> 6i/ e e£ 2 










61-62 












63-64 












65-66 



• CalcTBD : if this variable is set .TRUE, then in all chargino-, neutralino- and gluino 
decays the three body modes will be calculated. This option has to be taken with care, 
because it can slow down the program considerably. 

• ratioWoM : this variable is used to decide whether two body decays or three body decay 
modes are calculated in case of charginos, neutralino and gluino. The program tries 
first two-body decay modes. In the case that the ratio of the width T over the mass m 
of the decaying particle is small: T/m < ratioWoM, then three body decay modes are 
calculated. 

• Fgmsb and m32 : the F parameter and the gravitino mass in the GMSB model. These 
parameters are calculated from the input and are set to huge numbers in all other 
models. They are needed for the calculation of the decay width(s) of the NLSP into a 
gravitino. 

The next statements call the routine for the calculation of the cross sections provided 
L_CS = . TRUE . and kont = 0: 

If ( (L_CS) . and . (kont . eq . 0) ) then 
Call InitializeCrossSections (Ecms , Pm, Pp, ISR) 
Call CalculateCrossSections (Ecms , Pm, Pp, ISR & 
& , mSup, RSup, mf_u, mSdown, RSdown, mf_d, mglu & 
& , SigSup, SigSdown, mSlepton, RSlepton & 
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Table 4: Correspondence between the second indices for chargino partial widths (branching 
ratios) and the decay modes. 



mode 



index of gP_C ( BR_C) 



Xi 

Xi 

Xi 
Xi 
Xi 
Xi 

xi 
xi 
xi 
xi 



l m,k u m 

V ' l + 
^m h m 

dm,k 
^"m,k d m 

X°jW + 
X°H + 

xiz° 

xtA° 
Xth° 

xi H° 



1-6 
7-9 

10-15 
16-21 
22-25 
26-29 

30 

31 

32 

33 



Xi 
Xi 

X.i 

xi 
xi 
xi 
xi 



Xj ^"m d m 

Y°l + V 

A.j h m _rn 

9 ^"m d m 
Xl 

Xl dm dm 
Xl lm lm 
Xl ^m 



64-75 
76-87 
88-90 
91-93 
94-96 
97-99 
100-102 



& , mSneut, RSneut, SigSle, SigSn, mC, U, V, mN, N & 
& , SigC, SigChiO, mSO, RS0, vevSM, mPO, RPO, mSpm & 
& , RSpm, SigSO, SigSP, SigHp ) 

End If 

Here the following additional input is needed: 

• Ecms : the center of mass energy of the collider 

• Pm, Pp : degree of polarization of the incoming electron and positron, respectively 

• ISR : logical variable, if .TRUE, then initial state radiation is taken into account using 
the formulas given in [54] 

These variables can be set in the file CrossSections . in. The cross sections are stored in 
the variables starting with Sig which are summarized in Table El Please note, that in case 
of sfermions the structure of the variables is already put such that the case of generation 
mixing can easily be implemented. In the non-mixing case the cross sections are stored in 
the 2x2 diagonal blocks and they are sorted according to the generations as in the case of 
the sfermion mixing matrices. 

Finally the low energy constraints b — * 37, and Ap are calculated provided that 
calculation of the spectrum had been performed successfully (kont.eq.O): 

If (kont.eq.O) then 
Call CalculateLowEnergyConstraints (gauge , Y_l, Y_d, Y_u & 
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Table 5: Correspondence between the second indices for neutralino partial widths (branching 
ratios) and the decay modes. Note, that also charge conjugated states are given. 



mode 



index of gP_N ( BR_N) 



Y° 

At 

Y° 

At 

x\ 

Y° 

x\ 



r 1 

dm,k dm 

X°jA 

xjh° 



1-12 

13-18 

19-30 

31-42 

43-46 

47-50 
51-(24+i) 
(25+i)-(23+2 i) 
(26+i)-(22+3 i) 
(27+i)-(21+4 i) 



Y° 

Y° 

A.% 



Z°G 
h°G 



63 
64 
65 



Y° — 

At 

'X.i 

Y° — 

At 

Y° — 

At 

Y° — 

A^ 

Xi>l 
Xi>l 
Xi>l 

X?>i 

Xi>2 

Xi>2 
~0 

Xi>2 
Xi>2 

xl- 
xl- 
xl- 
xl- 



Xj Qm Qm 
Xj I'm 
G U m "Urn 
G dm dm 
~ * Xl U m 

— > y? d d 

Al u, m u, m 

-> Y? 1+ I 
Al L m L 

-> Xl^m 

* X2 ^"in V> m 
~^ X2 dm d m 
~* X.2 lm, lm 
~ * X2 

X3 ^m ^m 

X.3 d m d m 

Y° /+ /- 

A3 m "m 

X3 ^m 



rn 

Urn 



65 +j 
69 - 80 
81 - 92 
93 - 95 
96 - 98 
99 - 101 
102 - 104 
105 - 107 
108 - 110 
111 - 113 
114 - 116 
117 - 119 
120 - 122 
123 - 125 
126 - 128 
129 - 131 
132 - 134 
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Table 6: Correspondence between the indices for gluino partial widths (branching ratios) and 
the decay modes. Note, that also charge conjugated states are given. 



mode 



index of gP_G ( BR_G) 



9 
9 



dm,k d m 



1-12 
13-24 



9_ 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 



t± c 



25-26 



if 

Xl U m 

y° d d 

Al u m u m 
X2 U m u m 

X.2 dm d m 

X3 V>m 

X3 d m d m 

Xa dm d m 
Xi Qm Q m 

U W~ b 



26 +j 
31-33 
34-36 
37-39 
40-42 
43-45 
46-48 
49-51 
52-54 
55-66 
67-68 



& , mSpm2, RSpm, mC, U, V, mN, N , mSup2, RSup, mSdown2 & 
& , RSdown, mSlepton2, RSlepton, mSneut2, RSneut & 
& , BRbtosgamma, a_mu, Delta_Rho) 
Else 

BRbtosgamma = 0._dp 
a_mu = 0._dp 
Delta.Rho = 0._dp 
End If 

Here BRbtosgamma, a_mu, and Delta_Rho denote 10 4 x BR(b — > S7), the SUSY contributions 
to and the sfermion contributions to Ap, respectively. 
Afterward the statement 

Call WriteOutPutOCll, l.e-6_dp, l.e-3_dp) 

is used to write all information to the file connected with unit 11 (first entry). The second 
entry puts a lower bound on the branching ratios to be written. In the case above, branching 
ratios smaller than 10~ 6 will not be given. The third entry gives the minimum value for the 
cross section in fb which will be written to the output file. In the example given above cross 
sections smaller than 1CT 3 will not be written to the output file. 
The last statement closes all open files. 

call closingO ! closes the files 



8 Conclusions 

We have described SPheno, a program calculating the spectrum, branching ratios and cross 
sections of supersymmetric particle in e + e~ annihilation within the MSSM. The user can 
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Table 7: Correspondence between the indices for the partial widths (branching ratios) of the 
neutral Higgs bosons and the decay modes. The variables are gP_S0 (BR_S0) and gP_P0 
(BR_P0) for the partial decay widths (branching ratios) of the CP-even Higgs bosons (h°, 
H°) and CP-odd Higgs boson (A ). In case of gP_S0 ( BR_S0) the first (second) index denotes 
decay modes of h° (H°). Here stands for h°, H° and A . The index m runs from 1 to 3. 



mode 



h° 



A 









l + l~ 
i m l m 

dm dm 

Urn t^m. 



e 1 e 2 

— > e 2 62 

-+ fit fil 
fit fit 
-»• fit P>2 



'1 T 2 

-+ f 2 + 12 
— > U m U m 

-> di di 
dfdf 
-> d\ d 2 
-> §t sj 

S t S 2 
-> S 2 s 2 

bfbt 

-> Hf Ui 

-»• ct £r 

c 1 c 2 
-> c 2 c 2 

-> t 2 ^2 



1-3 
4-6 
7-9 



1-3 
4-6 
7-9 



1-3 
4-6 
7-9 



0- 

H° 
H° 

0- 

H° 

H° 

H° 
H° 
H° 

0" 

H° 
H° 

0" 
#° 

#° 

0- 

#° 

#° 

0" 
#° 

H° 

0- 
#° 

#° 

0- 

H° 



10 
11-12 

13 

14 
15-16 

17 

18 
19-20 

21 
21 + m 

25 
26-27 

28 

29 
30-31 

32 

33 
34-35 

36 

37 
38-39 

40 

41 
42-43 

44 

45 
46-47 

48 



11-12 



15-16 



19-20 



23-24 



27-28 



31-32 



35-36 



39-40 



43-44 



-»■ X° r X° s (r < s) 

<P^xtxi 

-»• xt X2 
0^X2X2 



49-58 

59 
60-61 

62 



49-58 

59 
60-61 

62 



H» - Z u Z y) 
H° -> W + W~ 
H° -> /i° /i° 
A -> /i° Z° 



46-55 

56 
57-58 

59 



63 
64 
70 



63 
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Table 8: Correspondence between the indices for the partial widths (branching ratios) of the 
charged Higgs and the decay modes. The variables are gP_Spm ( BR_Spm) and in case of 
partial decay widths (branching ratios). 



mode 



index 



H+ 



U 



m d m 



1-3 
4-6 



H+ 
H+ 
H+ 

H + 
H+ 

H + 
H+ 
H+ _> h° W~ 



el v e 

ffv T 
Ui dj 

Ci Sj 
Ar As 



7-8 
9-10 
11-12 

12 + 2*(j-l) + i 
16 + 2*0-1) + i 
20 + 2*(j-l) + i 
24 + 4*(r-l) + s 
34 



Table 9: Correspondence between the production cross sections and the variables used in the 
program. 



process 



Fortran name and type 



e 1 e 
e + e" 
e + e" 
e + e" 



Ui Uj 

di dj 
li lj 



real (dp) 
real (dp) 
real (dp) 
real (dp) 



SigSup(6,6) 
SigSdown(6,6) 
SigSle(6,6) 
SigSn(6,6) 



e 1 e 
e + e 



Afc An 



real (dp) 
real (dp) 



SigChiO(4,4) 
SigN(4,4) 



e 1 e 
e + e" 
e + e" 



H+H- 



real(dp) 
real (dp) 
real (dp) 



SigS0(2) 
SigSP(2) 
SigHp 
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choose between the following high scale models: minimal supergravity, minimal supergrav- 
ity including right handed neutrinos, gauge mediated supersymmetry breaking, anomaly 
mediated supersymmetry breaking, and string effective field theories based on 01 and Oil 
compactification. The calculation of the spectrum are done using two-loop renormalization 
group equations and the complete one-loop formulas for the SUSY masses. In case of the 
neutral Higgs bosons and the \i parameter leading two-loop effects are included. The masses 
and mixing angles are used to calculate the most important two body and three body decay 
modes. They are also used for the calculation of the SUSY production cross sections in e + e~ 
annihilation. Here the effects of initial state radiation and longitudinal beam polarization 
is included. Finally the following low energy quantities are calculated: BR(b — > sj), the 
super symmetric contributions to the anomalous magnetic moment of the muon and the 
sfermion contributions to the p parameter. Starting with version 2.2.0 SPheno allows for 
input and output according to the SUSY Les Houches accord [2~5| . 

The program is set up in such a way that extensions can easily be included. The plans 
for upcoming versions are to include complex phases for the supersymmetric parameters, 
to include generation mixing, to include QCD and Yukawa corrections for various processes 
such as Sfermion and Higgs production and decays. In addition beam strahlung for various 
collider designs will be implemented. 
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A Switches 

In this appendix we describe the switches for influencing the behaviour of the program. Inside 
the main program one can set two scales: 

1. The electroweak scale M EWS b, which is the scale where the loop contributions to the 
masses and mixing matrices are calculated. The default is to calculate this scale from 
M EWS b = y/m^m^. By calling 

Call SetRGEScale(l.e3_dp**2) 

Mew sb will be set to the fixed value of 10 3 GeV in this example. Note that the input 
is the scale squared. In the case one uses a zero or a negative number as input for 
SetRGEScale then the scale will be calculated from the stop masses. 
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2. The high energy scale, where the boundary conditions of the model under study are 
set. Except for GMSB, where the scale is fixed by default, the high scale is calculated 
from the requirement g\ — g^. By calling 

Call SetGUTScale(2.el6_dp) 

Mqut w iU be set to 2 ■ 10 16 GeV in this example. In the case one uses a zero or a 
negative number as input for SetGUTScale then the scale will be calculated from the 
requirement g\ = g 2 except for GMSB. 

In general the strong coupling g s will be different from gi and g 2 in GUT theories if one works 
at the two loop level [ZZj. In case someone wants to enforce strict universality at the high 
scale, this can be done by using the following statement: 

test = SetStrictUnification(.TRUE.) 

where test is a logical variable to which the old value of SetStrictUnification is assigned. 
The same functions can also be used to set this option . FALSE . . 

Starting with version 2.2.0 it also possible to set these values using the SUSY Les Houches 
accord as described in Appendix O 

B Input files 

In this section the input files are described. Among these files only the file HighScale . in has 
to be provided by the user. The other files Control . in, CrossSections . in and StandardModel . in 
can be used to change the default values which are given below in Appendices IB. 11 IB. 21 and 
IB. 41 An alternative way to provide SPheno with input is by the SUSY Les Houches accord 
as described in Appendix O 

B.l Control. in 

This file contains three entries as shown below: 
! ErrorLevel 

.True. ! Calculation of branching ratios 

.True. ! Calculation of cross sections 

The values given above are the default values inside the program and are used if the the 
file Control, in is not present. Here ErrorLevel is an integer in the range [-2,2] where the 
numbers correspond to the following behaviour of the program: 

-2 do not print severe warnings 

-1 do not print warnings 

print every warning 

1 abort in case of a severe warning 

2 abort even in case of a warning 

A warning is called severe if either a result is unphysical or if a numerical procedure gives an 
unreliable result. 
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B.2 CrossSections.in 



This file contains four entries as shown below: 



500. 

0. 

0. 

.True 



degree of longitudinal polarization of electrons 
degree of longitudinal polarization of positrons 
calculation of initial state radiation if .TRUE. 



c.m.s. energy in GeV 



The values given above are the default values inside the program and are used if the the file 
CrossSections.in is not present. 



In this section we describe the input file for the high scale boundary conditions. The package 
contains several files starting with Highscale . in and ending in the models described below. 
One has to rename the model file to the name Highscale. in to use it as input for SPheno. 
Note that in all examples below the value of Aq given below will be multiplied by the Yukawas 
coupling at the high scale and that this product enters the RGEs. For most of the examples 
below we have used the so-called SPS points defined in [75] . In Appendix |D] we display the 
output for the point SPSla. 

B.3.1 mSUGRA 

The minimal SUGRA version is defined by four parameters and the sign of the \x parameter. 
The parameters are the gaugino mass parameter Mi/ 2 , the scalar mass parameter M , the 
trilinear parameter A as well as tan /3. The file reads for example 



B.3.2 mSUGRA including right handed neutrinos 

In this case one needs four more input values compared to the case of mSUGRA described 
above: a common right handed neutrino mass m UR and the light neutrino masses m Vi (i = 
1, 2, 3). The file reads for example 



B.3 HighScale.in 



mSugra 
250. 
100. 

-100. 
10. 



M.l/2 

M_0 

A_0 

tan (beta) 
sign of mu 

if 2-loop RGEs should be used 



1. 
. TRUE . 



mSugra 
250. 
100. 

-100. 



M.l/2 

M_0 

A_0 

tan (beta) 
sign of mu 
m_nu_R 



10. 
1. 

I.el4 

l.e-14 3.e-12 0.06e-9 
. TRUE . 



m_nu_i 



if 2-loop RGEs should be used 



30 



B.3.3 GMSB 



The GMSB scenario is characterized by two mass parameters Mm and A; the multiplicity 
N 5 and N w of messengers in the 5 + 5 and 10 + 10 vector-like multiplets, respectively; tan /3 
and the sign of \i as described in Sect. 12.3.31 In addition one can set a common value for 
the A parameters at the scale Mm- Note, however, that in the minimal model this value is 
practically zero. The file reads for example 

GMSB 
100000. 
200000. 
1 

0. 
15. 
1. 

. TRUE . 



Lambda 

M_M 

N_5 

N_10 

AO 

tan (beta) 
sign of mu 

if 2-Iood RGEs should be used 



B.3.4 AMSB 

The implemented AMSB scenario is characterized by the gravitino mass 777,3/2, a common 
scalar mass Mq, tan/3 and the sign of \x. The file reads for example 



AMSB ! model 

60000. ! M.3/2 

450. ! M_0 

10. ! tan(beta) 

1 . ! sign of mu 

.TRUE. ! if 2-loop RGEs should be used 



B.3.5 String I 

This scenario is characterized by the gravitino mass m 3 / 2 , the common vev < t > of the 
moduli fields, the string coupling squared g 2 s1 the sine squared of the mixing angle between the 
dilaton fields and moduli fields sin 2 9, the parameter 5gs of the Green-Schwarz counter-term, 
the modular weights rii characterizing the couplings between moduli fields and matter fields, 
which are assumed to be generation independent in the current implementation. Moreover, 
one needs to specify tan (3 and the sign of \l. The file reads for example as 



String.OI 




180. 


! M_3/2 


14.0 


! <t> 


0.5 


! g-s~2 


0.9 


! sin~2(theta) 


0. 


! delta.GS 


-1 -3 


! n_E n_L 


1-2 


! n_D n_U n_Q 


-1 -1 


! n_Hl n_H2 


10. 


! tan(beta) 
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-1. 

. TRUE . 



! phase (mu) 

! if 2-loop RGEs should be used 



B.3.6 String II 

This scenario is characterized by the gravitino mass m.3/2, the common vev of the moduli 
fields < t >, the string coupling squared g%, the sine squared of the mixing angle between 
the dilaton fields and moduli fields sin 2 9, the parameter of the Green-Schwarz counter- 
term 5gs- Moreover, one needs to specify tan/3 and the sign of fi. There are two different 
scenarios implemented denoted as String_OIIa and String_OIIb corresponding to boundary 
conditions (A) and (B) of [HI], respectively. The files read for example as 

String.OIIa 



300. ! M.3/2 

14.6 ! <t> 

0.5 ! g_s~2 

0.9 ! sin~2(theta) 

0. ! delta.GS 
5. ! tan (beta) 

1. ! phase (mu) 

.TRUE. ! if 2-loop RGEs should be used 

and 

String.OIIb 

300. ! M.3/2 

14.6 ! <t> 

0.5 ! g_s~2 

0.9 ! sin~2(theta) 

0. ! delta.GS 
5. ! tan (beta) 

1. ! phase (mu) 

.TRUE. ! if 2-loop RGEs should be used 



B.3.7 SUGRA 

This input files serves as interface for more general models with gauge couplings unifica- 
tion. Here the user can specify non-universal gaugino masses at the high scale Mi/ 2 [U(1)], 
Mi/2 [577(2)], Mi/ 2 [SU(3)], 15 different values of the sfermion mass parameters for every type 
of sfermions: M% .., M~ .., M% .., M% .., M% ..; two Higgs mass parameters M% and M% ; 
nine different A parameters A ^ tii , A ^,u and A ^u. Here ii denotes that only the diagonal 



entries can be set. 




Sugra 






480. 300. 


300. 


! M_l/2_i 


150. 150. 


150. 


! M0_E_ii 


150. 150. 


150. 


! M0_L_ii 


150. 150. 


150. 


! M0_D_ii 


150. 150. 


150. 


! M0_Q_ii 



32 



150. 150. 
150. 150. 
0. 0. 0. 
0. 0. 0. 
0. 0. 0. 
10. 
1. 

. TRUE . 



150. 



M0_U_ii 

M0_H_i 

A0_u_ii 

A0_d_ii 

A0_e_ii 

tan (beta) 

phase (mu) 

if 2-loop RGEs should be used 



B.3.8 MSSM 

This input file serves for the cases that one wants to start with low energy parameters for the 
calculation of masses, decays and/or production rates. The input consists of the electroweak 
gaugino mass parameters M\ and M2, the gluino mass rrig, 15 sfermion mass parameters 

Aan; tan/5. Moreover, 



Mjz M , M Ltii , 



M D,iv M u,iv M Q,iv 9 trilinear parameters 



one has to specify the renormalisation scale Q where all the parameters are given as well as 
fj, and the mass of the pseudoscalar Higgs m^. 



MSSM 

99.13 192.74 580.51 
136.23 136.23 133.55 
196.64 196.64 195.75 
519.53 519.53 516.86 
539.86 539.86 495.91 
521.66 521.66 424.83 
0. 0. -510.01 
0. 0. -772.66 
0. 0. -254.20 
10. 454.65 
352.39 393.63 



M_l M_2 m_gluino 

M_E_i 

M_L_i 

M_D_i 

M_Q_i 

M_U_i 

A_u 

A_d 

A_e 

tan (beta) Q 
mu m A 



In the case that the keyword MSSMtree is used instead of MSSM the masses are calculated 
using tree-level formulas instead of loop corrected masses. The renormalisation scale Q is 
absent in this case, e.g. the two last line read as 



10. 

352.39 



393.63 



! tan (beta) 
! mu m_A 



B.4 StandardModel.in 

This file contains the values of the Standard Model parameters and must include all lines 
given below. Otherwise the default values given in the listing below are used: 



91.1876 
2.4952 

0.0336 0.0336 0.037 

0.2 

2.118 

0.1 0.1 0.1 



m_Z 

Gamma.Z 

Br(Z -> 1 1) 

Br(Z -> invisible) 

Gamma.W 

Br(W -> 1 nu) 
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0.51099890e-3 

0.105658357 

1.7770 

2.00 

0.003 

1.2 

175.0 

0.007 

0.12 

4.25 

137.0359998 

0.1172 

1.16639e-5 

0.224 

0.0413 

0.00363 

0. 

2.19709e-6 
3.4e-13 



scale Q where the masses of the light quarks u,d,s,c are g 

m_u(Q) 

m_c(Q) 

m_t, pole mass 

m_d(Q) 

m_s(Q) 

m_b(m_b) 

1 . / alpha 

alpha_s(m_Z) 

G_F, Fermi constant 
sl2 of CKM particle data book 1998, 90%, 0.217-0.222 
s23 of CKM 0.036-0.042 
! sl3 of CKM 0.0018 - 0.0044 
phase of CKM 0.-2 Pi 



life time of muon 
life time of tau 



m_e 



m_mu 



m_tau 



All masses are given in GeV and Br denotes "branching ratio" in the list above. 

C Implementation of SUSY Les Houches Accord 

Starting with version 2.2.0 SPheno allows for input and output according to the SUSY Les 
Houches accord [2H]. The name of the input file is LesHouches . in and the output will be 
written to the file SPheno . spc. In the following we summarize unsupported features as well 
extensions of this standard. The unsupported features are: 

• a\ m {rnz) MS of the block SMINPUTS is ignored as input. Instead the value a em (0) can 
be changed in the block SPhenoInput described below. The corresponding value for 
the output will be given. 

• In block EXTPAR the entries 51-53 are ignored as the corresponding formulas are not 
(yet) implemented in SPheno. 

• Currently there is no information concerning warnings and errors in block SPINF0. This 
will be changed within the one of the next versions. Please check the file Messages . out 
for this information. 

The current implementation requires that the block M0DSEL is read in before the block EXTPAR 
is read in. In the case that an unknown entry appears, a warning message is printed and 
SPheno tries to proceed. SPheno stops execution in the case that the model input is not 
complete. 

For the SPheno specific input the block SPhenoInput has to be used. In this block 
switches can be set, SM input beside the one of the block SMINPUTS can be set. Moreover, 
the information for the cross section calculation can be given here. Within this block the 



34 



following flags and parameters can be set, with general structure id value: 

1 : setting the error level as described Appendix IB. II 

11 : if value=l (0) then (no) branching ratios are calculated 

12 : only branching ratios larger than value are written out 

21 : if value=l (0) then (no) cross sections are calculated 

22 : cms energy for e + e~ annihilation 

23 : value gives degree of polarisation for e~ beam 

24 : value gives degree of polarisation for e + beam 

25 : only cross sections larger than value (in fb) are written out 

31 : a fixed value for the GUT scale is used if value is larger than 

32 : if value=0 then g3(m GUT ) can be different from gi{m GUT ) = ^(^gi/t); 

if value=l then strict unification giimaur) — 92( m GUT) — 9z{. m GUT) is enforced 

33 : a fixed value for the renormalization scale Qewsb is used if value is larger than 

41 : sets value of Z-boson width 

42 : sets value of W-boson width 

51 : sets value of electron mass 

52 : sets value of muon mass 

61 : sets scale where the running masses for light quarks (u, d, s, c) are defined 

62 : sets value of u-quark mass 

63 : sets value of c-quark mass 

64 : sets value of d-quark mass 

65 : sets value of s-quark mass 

Here is an example for this block with the default values of SPheno 



Block SPhenoInput 


# 


SPheno specific input 


1 


-1 




# 


error level 


11 


1 




# 


calculate branching ratios 


12 


1 


00000000E-04 


# 


write only branching ratios larger than this value 


21 


1 




# 


calculate cross section 


22 


5 


00000000E+02 


# 


cms energy in GeV 


23 





00000000E+00 


# 


polarisation of incoming e- beam 


24 





00000000E+00 


# 


polarisation of incoming e+ beam 


25 


1 


00000000E-04 


# 


write only cross sections larger than this value [fb] 


31 


-1 


00000000E+00 


# 


m_GUT, if < than it determined via g_l=g_2 


32 







# 


require strict unification g_l=g_2=g_3 if ' 1' is set 


33 


-1 


00000000E+00 


# 


Q.EWSB, if < than Q_EWSB=sqrt (m_~tl m_~t2) 


41 


2 


49520000E+00 


# 


width of the Z-boson 


42 


2 


11800000E+00 


# 


width of the W-boson 


51 


5 


10998900E-04 


# 


electron mass 


52 


1 


05658357E-01 


# 


muon mass 


61 


2 


00000000E+00 


# 


scale where quark masses of first 2 gen. are defined 


62 


3 


00000000E-03 


# 


m_u(Q) 


63 


1 


20000000E+00 


# 


m_c(Q) 


64 


7 


00000000E-03 


# 


m_d(Q) 


65 


1 


20000000E-01 


# 


m_s(Q) 



For the output the extensions below have been defined. The information concerning the 
cross section is written out using a SPheno specific block called SPhenoCrossSections. The 
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first line of this block gives the information on the cms energy, the polarization of the incoming 
beams as well if ISR is included or not, for example for = 500 GeV and unpolarized beams: 

Block SPhenoCrossSections # cross sections 

XS 11 -11 500.0 0.00 0.00 1 # e+ e- XS, Pe-, Pe+, including ISR 

The FORTRAN format is in this case: 

Format ("XS 11 -11 ",F7.1," " ,F5.2," \F5.2," " , A) 

The cross sections (in fb) themself are written as 

# Sigma [fb] NDA ID1 ID2 

2.83574498E+02 2 2000011 -2000011 # ~e_R- ~e_R+ 

Here the first entry gives the cross section in fb, the second entry specifies the number of 
produced particles, the subsequent two integers give the PDG code of the particles. We have 
used the FORTRAN format 

Format(3x,lP,el6.8,0p,3x,I2,3x,2(i9,lx) ,2x," # ",A) 

As an example we give the cross sections for the SPSla scenario at an 500 GeV e + e~ linear 
collider with unpolarized beams: 

Block SPhenoCrossSections # cross sections 

XS 11 -11 500.0 0.00 0.00 1 # e+ e- XS, Pe-, Pe+, including ISR 





Sigma [fb] 


NDA 


ID1 


ID2 








2 


.83574498E+02 


2 


2000011 


-2000011 


# 


~e_R- 


~e_R+ 


7 


.79728001E+01 


2 


2000011 


-1000011 


# 


~e_R- 


~e_L+ 


4 


.57495061E+01 


2 


1000011 


-1000011 


# 


~e_L- 


~e_L+ 


5 


.47916441E+01 


2 


2000013 


-2000013 


# 


~mu_R- 


~mu_R+ 


6 


.00045490E-03 


2 


2000013 


-1000013 


# 


~mu_R- 


~mu_L+ 


1. 


.90114309E+01 


2 


1000013 


-1000013 


# 


~mu_L- 


~mu_L+ 


5 


.96228076E+01 


2 


1000015 


-1000015 


# 


~tau_l- 


~tau_l+ 


1. 


. 26426385E+00 


2 


1000015 


-2000015 


# 


~tau_l- 


~tau_2+ 


1. 


.59684572E+01 


2 


2000015 


-2000015 


# 


~tau_2- 


~tau_2+ 


4 


.52889205E+02 


2 


1000012 


-1000012 


# 


~nu_eL 


~nu_eL* 


1. 


.36168303E+01 


2 


1000014 


-1000014 


# 


~nu_muL 


~nu_muL* 


1 


.39168830E+01 


2 


1000016 


-1000016 


# 


~nu_tauL 


~nu_tauL* 


2 


.75869582E+02 


2 


1000022 


1000022 


# 


chi_10 chi. 


.10 


6 


.56937491E+01 


2 


1000022 


1000023 


# 


chi_10 chi. 


.20 


7 


. 10141133E+00 


2 


1000022 


1000025 


# 


chi_10 chi. 


.30 


8 


.27993814E-01 


2 


1000022 


1000035 


# 


chi_10 chi. 


.40 


6 


.90281358E+01 


2 


1000023 


1000023 


# 


chi_20 chi. 


.20 


1 


. 60903760E+02 


2 


1000024 


-1000024 


# 


chi_l- chi. 


.1+ 


2 


.47077869E+01 


2 


25 


23 


# 


hO Z 





The information concerning the value of low energy observables (BR(b — > s'y), SUSY 
contribution to (g — 2)^ and A(p)) is written out using a SPheno specific block called 
SPhenoLowEnergy. We use the following identifier: 
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1 : BR(6 -> 57) 

2 : SUSY contribtutions to (g — 2) M 

3 : SUSY contribtutions to A(p) 

As an example we give here the output for the SPSla scenario: 

Block SPhenoLowEnergy # low energy observables 

1 4.55809155E+00 # BR(b -> s gamma) 

2 5.42193822E-09 # (g-2)_muon 

3 1.97608480E-04 # Delta(rho) 



D Sample output 

Here we give the content of the file SPheno . out provided one uses the content of HighScale . in 
for the mSUGRA scenario described in Appendix IB. 3. II and the default values of the files 
Control . in, CrossSections . in and StandardModel . in. 

SPheno output file 

Version 2.2.0 , created: 11.03.2004, 21:06 



mSugra input at the GUT scale 2 . 3626654756304632E+16 

M.l/2 : 2.5000000000000000E+02 

M_0 : 1.0000000000000000E+02 

A_0 : -1.0000000000000000E+02 

tan(beta) at m_Z : 10.0000000000000000 

phase (mu) : 1.0000000000000000 



Parameters at the scale 



4 . 8391963443672887E+02 



g' S g- 3 

3.61193779E-01 6.46479958E-01 1 . 09493111E+00 

Y_e Y_mu Y_tau 

2.88471294E-05 5 . 96467232E-03 1 . 00322460E-01 



Y_u Y_c Y_t 

8.84695152E-06 3 . 53878053E-03 8 . 92107505E-01 

Y_d Y_s Y_b 

1.91375025E-04 3 . 28071453E-03 1 . 37357260E-01 

Gaugino mass parameters 

1 . 0165957170339030E+02 1 . 9176076501787236E+02 5 . 8493277901867339E+02 

mu, B 

3 . 5735299500991079E+02 1 . 6699996261573800E+04 



Slepton mass parameters 
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A_l 

-2 . 5325613989602988E+02 
M2_E 

1 . 844076826831 1054E+04 
M2_L 

3 . 8191007905480838E+04 

Squark mass parameters 
A_d 

-8 . 5467574322023552E+02 
A_u 

-6 . 7949078243535337E+02 
M2_D 

2 . 744772242274885 1E+05 
M2_U 

2 . 7672441991504660E+05 
M2_Q 

2 . 9674725559000188E+05 



-2 . 5325004783357019E+02 
1 . 8438451866710398E+04 
3 . 8189867711691244E+04 



-2 . 5153292012959056E+02 
1 . 7787517468332015E+04 
3 . 7869523746231192E+04 



-8 . 5467228770510383E+02 -7 . 9058921072895021E+02 

-6 . 7948707229660602E+02 -4 . 9657599289477690E+02 

2 . 7447520834575844E+05 2 . 71 14817843756522E+05 

2 . 7672246826518845E+05 1 . 7613291366656474E+05 

2 . 9674529488349473E+05 2 . 4602686847035179E+05 



Higgs mass parameters 

3 . 2530345774801379E+04 -1 . 2784263452390164E+05 



Masses and mixing matrices 

Gluino : 6 . 0402198802780606E+02 



1.0000000000000000 



Charginos 

1 . 7973283139106542E+02 
U 

-0.91793 0.39673 
0.39673 0.91793 

V 

-0.97122 
0.23819 



3 . 8228031192319548E+02 



0.23819 
0.97122 



Neutralinos 
97.1655612758851674 
N 



1 . 8073254909635224E+02 3 . 6466658683987458E+02 3 .81926175864457491 



1 


1 ( 


-0.98576, 


0.00000) 


1 


2 ( 


0.05614, 


0.00000) 


1 


3 ( 


-0.14885, 


0.00000) 


1 


4 ( 


0.05449, 


0.00000) 


2 


1 ( 


-0.10391, 


0.00000) 


2 


2 ( 


-0.94275, 


0.00000) 


2 


3 ( 


0.27483, 


0.00000) 


2 


4 ( 


-0.15776, 


0.00000) 


3 


1 ( 


0.00000, 


0.06048) 


3 


2 ( 


-0.00000, 


-0.09026) 
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3 3 ( 

3 4 ( 

4 1 ( 
4 2 ( 
4 3 ( 
4 4 ( 



-0.00000, 
-0.00000, 

0.11757, 
-0.31610, 
-0.64770, 

0.68319, 



e-sneutrino mass 
mu-sneutrino mass 
tau-sneutrino mass 



-0.69484) 
-0.71091) 
0.00000) 
0.00000) 
0.00000) 
0.00000) 

1 . 9127930364369408E+02 
1 . 9127614461238977E+02 
1 . 9038540593550744E+02 



selectron masses 

1 . 4393792438692086E+02 2 . 0712432781915908E+02 
R_e 

0.00009 1.00000 
-1.00000 0.00009 

smuon masses 

1 . 4390451752393614E+02 2 . 0713856716885780E+02 
R_mu 

0.01790 0.99984 
-0.99984 0.01790 

stau masses 

1 . 3484343437971901E+02 2 . 1068049221898704E+02 
R_tau 

0.26425 0.96445 
-0.96445 0.26425 

u-squark masses 

5 . 4756917879839784E+02 5 . 6524131281347229E+02 
R_u 

0.00006 1.00000 
-1.00000 0.00006 

c-squark masses 

5 . 4755738858800510E+02 5 . 6524965433179352E+02 
R_c 

0.02389 0.99971 
-0.99971 0.02389 

t-squark masses 

3 . 9945389533787068E+02 5 . 8628028438847809E+02 
R_t 

0.55340 0.83291 
-0.83291 0.55340 
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d-squark masses 

5 . 4730506514994352E+02 5 . 7064687850962798E+02 
R_d 

0.00058 1.00000 
-1.00000 0.00058 

s-squark masses 

5 . 4730068877040878E+02 5 . 7064721260745614E+02 
R_s 

0.01000 0.99995 
-0.99995 0.01000 

b-squark masses 

5 . 1512514815767906E+02 5 . 4710587491028832E+02 
R_b 

0.94697 0.32132 
-0.32132 0.94697 

m_A0, m_H+ 

3 . 9942973783298265E+02 4 . 0772744504014719E+02 
m_h0, m_H0 

1 . 1080888761 166447E+02 3 . 9980703384665037E+02 
R_S0 

0.11371 0.99351 
-0.99351 0.11371 



Low energy constraints 
10~4 Br(b -> s gamma) 
Delta(a_mu) 
Delta(rho) 



0.4558514E+01 
0.5415330E-08 
0.1959489E-03 



Anti particles are marked with a * in case of 
(s)neutrinos and (s) quarks in the decay section. 

Decay widths (GeV) and branching rat 

Selectron_l 

Neutralino.l e 0.21510088 100.00000000 

Total width : 0.21510088 



Selectron_2 
Neutralino_l e 
Neutralino_2 e 
Chargino_l neutrino 



0. 12813720 
0.04910231 
0.08853247 



48.21320960 
18.47535119 
33.31143921 
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Total width 



0.26577199 



Smuon_l 

Neutralino.l mu 0.21484110 100.00000000 

Total width : 0.21484110 



Smuon_2 

Neutralino. 

Neutralino. 

Chargino_l 
Total width 



1 mu 

2 mu 
neutrino 



0. 12829264 
0.04912578 
0.08857176 
0.26599018 



48.23209605 
18.46902142 
33.29888253 



Stau.l 

Neutralino.l tau 0.15075247 100.00000000 

Total width : 0.15075247 



Stau_2 

Neutralino. 

Neutralino. 

Chargino_l 
Total width 



1 tau 

2 tau 
neutrino 



. 16323638 
0.05485389 
0.09802725 
0.31611751 



51.63787855 
17.35237248 
31.00974897 



e-Sneutrino 
Neutralino. 
Neutralino. 
Chargino_l 

Total width 



1 neutrino 

2 neutrino 
e 



0. 16123722 
0.00715724 
0.02056576 
0.18896022 



85.32866208 
3.78769424 
10.88364368 



mu-Sneutrino 
Neutralino.l neutrino 
Neutralino_2 neutrino 
Chargino_l mu 

Total width : 



0.16123086 85.33568036 

0.00715318 3.78600969 

0.02055318 10.87830995 
0.18893721 



tau-Sneutrino 
Neutralino.l neutrino 
Neutralino_2 neutrino 
Chargino_l tau 

Total width : 



0.15943476 87.31209427 

0.00605108 3.31378333 

0.01711746 9.37412241 
0.18260329 
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Sdown.l 

Neutralino. 

Neutralino. 

Neutralino. 

Neutralino. 

Chargino_l 
Total width 



1 d-quark 

2 d-quark 

3 d-quark 

4 d-quark 
u- quark 



0.28770802 
0.00270653 
0.00035708 
0.00114841 
0.00000154 
0.29192159 



98.55660629 
0.92714381 
0.12232218 
0.39339795 
0.00052910 



Sdown_2 

Neutralino. 

Neutralino. 

Neutralino. 

Neutralino. 

Chargino_l 

Chargino_2 
Total width 



1 d-quark 

2 d-quark 

3 d-quark 

4 d-quark 
u- quark 
u- quark 



0.12854555 
1 . 63736438 
0.00855909 
0.08259529 
3.24403336 
0.22691457 
5.32801224 



2.41263613 
30.73124285 

0.16064329 

1.55020837 
60.88637215 

4.25889723 



S-strange_l 
Neutralino. 
Neutralino. 
Neutralino. 
Neutralino. 
Chargino_l 
Chargino_2 

Total width 



1 s -quark 

2 s-quark 

3 s-quark 

4 s-quark 
c- quark 
c- quark 



0.28771851 
0.00292980 
0.00036645 
0.00114683 
. 00045378 
0.00000057 
0.29261594 



98.32632702 
1.00124494 
0.12523375 
0.39192166 
0.15507723 
0.00019540 



S-strange_2 
Neutralino. 
Neutralino. 
Neutralino. 
Neutralino. 
Chargino_l 
Chargino_2 

Total width 



1 s-quark 

2 s-quark 

3 s-quark 

4 s-quark 
c- quark 
c- quark 



0.12853671 
1.63714210 
0.00858872 
0.08262979 
3.24352858 
0.22701084 
5 . 32743673 



2.41273081 
30.73039024 

0.16121673 

1.55102333 
60.88347439 

4.26116451 



Sbottom.l 
Neutralino_l 
Neutralino_2 
Neutralino_3 
Neutralino_4 
Chargino_l t- 
Stop.l W- 



b-quark 
b-quark 
b-quark 
b-quark 
quark 



0.16674913 
1 . 34234704 
0.01963712 
0.04237541 
1.72934949 
0.56211179 



4.31705130 
34.75269181 
0.50839523 
1.09707817 
44.77199120 
14.55279229 
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Total width 






3. 


.86256998 






Sbottom_2 














Neutralino. 


.1 


b-quark 


0. 


.23969032 


31. 


.70912297 


Neutralino. 


.2 


b-quark 


0. 


.09310073 


12. 


.31648674 


Neutralino. 


.3 


b-quark 


0, 


.04221537 


5, 


.58475713 


Neutralino. 


.4 


b-quark 


0. 


.05849015 


7, 


.73778189 


Chargino_l 


t- 


-quark 


0. 


. 12448615 


16. 


.46852800 


Stop.l W- 






0, 


.19792061 


26. 


18332328 


Total width 






0. 


.75590333 






Sup_l 














Neutralino. 


.1 


u-quark 


1. 


.15145952 


98. 


.55637070 


Neutralino. 


.2 


u-quark 


0. 


.01083087 


0. 


.92704170 


Neutralino. 


.3 


u-quark 


0. 


.00143110 


0. 


. 12249171 


Neutralino. 


.4 


u-quark 


0, 


.00460431 


0. 


. 39409473 


Total width 






1. 


. 16832581 






Sup_2 














Neutralino. 


.1 


u-quark 


0. 


.03594243 


0. 


.65176659 


Neutralino. 


_2 


u-quark 


1. 


.75313101 


31. 


.79062555 


Neutralino. 


.3 


u-quark 


0. 


.00499759 


0. 


.09062447 


Neutralino. 


.4 


u-quark 


0. 


.06006824 


1. 


.08925513 


Chargino_l 


d- 


-quark 


3 


.58197347 


64. 


.95417438 


Chargino_2 


d- 


-quark 


0, 


.07850353 


1. 


.42355387 


Total width 






5. 


.51461627 






S-charm_l 














Neutralino. 


.1 


c-quark 


1. 


.15081537 


98. 


.28078206 


Neutralino. 


.2 


c-quark 


0, 


.01188610 


1. 


.01508491 


Neutralino. 


.3 


c-quark 


0. 


.00147943 


0. 


. 12634444 


Neutralino. 


.4 


c-quark 


0, 


.00458850 


0. 


.39186267 


Chargino_l 


s- 


-quark 


0. 


.00217691 


0. 


.18591022 


Total width 


• 




1. 


. 17094649 






S-charm_2 














Neutralino. 


.1 


c-quark 


0. 


.03656941 


0. 


.66343674 


Neutralino. 


.2 


c-quark 


1. 


.75204300 


31. 


.78529687 


Neutralino. 


.3 


c-quark 


0. 


.00499658 


0. 


.09064717 


Neutralino. 


.4 


c-quark 


0. 


.06014129 


1. 


.09107406 


Chargino_l 


s- 


-quark 


3. 


.57978849 


64. 


.94397663 


Chargino_2 


s- 


-quark 


0. 


.07857902 


1. 


.42556852 
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Total width 



5.51211779 



Stop.l 
Neutralino_l t-quark 
Neutralino_2 t-quark 
Chargino_l b-quark 
Chargino_2 b-quark 
c-quark neutralino_l 
c-quark neutralino_2 
W b-quark neutralino_l 

Total width : 



0.39740101 
0.24141364 
1.39258046 
0.02117232 
0.00039523 
0.01609207 
0.00000851 
2.06906324 



19.20680823 
11.66777483 
67.30487641 
1.02328034 
0.01910206 
0.77774681 
0.00041131 



Stop_2 
Neutral ino. 
Neutralino. 
Neutralino. 
Neutralino. 
Chargino_l 
Chargino_2 
Stop.l Z 
Stop.l hO 

Total width 



1 t-quark 

2 t-quark 

3 t-quark 

4 t-quark 
b-quark 
b-quark 



0.22065736 
0.63812557 
0.31053610 
1.44635779 
1.59111038 
1 . 47358454 
1.38482788 
0.26895379 
7.33415341 



3.00862752 
8.70073930 
4.23410966 
19.72085545 
21.69453369 
20.09208783 
18.88190500 
3.66714154 



Chargino_l 
Smuon_l neutrino 
Stau_l neutrino 
Neutralino_l W 
neutralino_l d~* u 
neutral ino_l s~* c 
neutralino_l e~+ nu 
neutralino_l mu~+ nu 
neutralino_l tau~+ nu 

Total width : 



Chargino_2 
Selectron_2 neutrino 
Smuon_2 neutrino 
Stau_l neutrino 
Stau_2 neutrino 
e-sneutrino e 
mu-sneutrino mu 
tau-sneutrino tau 
Neutralino.l W 
Neutralino 2 W 



0.00003886 
0.01301743 
0.00055115 
0.00000003 
0.00000003 
0.00002675 
0.00002675 
0.00002661 
0.01368761 



0.28387005 
95.10376372 
4.02662386 
0.00020879 
0.00020863 
0.19544840 
0.19544086 
0.19442906 



0.12483567 
0. 12487681 
0.00129050 
0.13608583 
0.05066962 
0.05073214 
0.06846536 
0. 17054406 
0.71196465 



5.01198171 
5.01363320 
0.05181199 
5.46366024 
2.03431606 
2.03682600 
2 . 74879077 
6.84711094 
28.58440836 
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Chargino_l Z 
Chargino_l hO 
neutralino_l b~* 
neutralino_2 d~* 
neutralino_2 s~* 
neutralino_2 b~* 
chargino_l u u~* 
chargino_l c c~* 
chargino_l d d~* 
chargino_l s s~* 
chargino_l b b~* 
Total width : 



0.60191556 
0.44847615 
0.00028615 
0.00000457 
0.00000457 
0.00002663 
0.00001046 
0.00001045 
0.00000445 
0.00000446 
0.00052754 
2 . 49074476 



24. 16608771 
18.00570494 
0.01148873 
0.00018361 
0.00018367 
0.00106920 
0.00041986 
0.00041949 
0.00017860 
0.00017891 
0.02117995 



Neutralino_l : stable 



Neutralino_2 



Selectron~-_: 


1 e~+ 




0. 


.00067750 


3. 


.46314577 


Selectron~+_: 


1 e~- 




0. 


.00067750 


3. 


.46314577 


Smuon~-_l mu 


~+ 




0. 


.00070088 


3. 


.58266441 


Smuon~+_l mu 






0. 


.00070088 


3. 


.58266441 


Stau~-_1 tau 


~+ 




0. 


.00836773 


42. 


.77287626 


Stau"+_1 tau 






0. 


.00836773 


42. 


.77287626 


Neutralino_l 


photon 




0. 


.00000004 


0. 


.00019852 


neutral ino_l 


u u~* 




0. 


.00000350 


0. 


.01786695 


neutral ino_l 


c c~* 




0. 


.00000349 


0. 


.01784001 


neutral ino_l 


d d~* 




0. 


.00000463 


0. 


.02368626 


neutral ino_l 


s s~* 




0. 


.00000463 


0. 


.02368617 


neutral ino_l 


b b~* 




0, 


.00000476 


0. 


. 02435045 


neutral ino_l 


nu_e 


nu_e~ 


0. 


.00001156 


0. 


.05908396 


neutral ino_l 


nu_mu 


nu_mu 


0. 


.00001156 


0. 


.05909905 


neutral ino_l 


nu_tau 


nu_ta 


0. 


.00001244 


0. 


.06357803 


neutral ino_l 


e"- 


e~+ 


0, 


.00000506 


0. 


.02584007 


neutral ino_l 


mu~- 


mu~+ 


0. 


.00000505 


0. 


.02582376 


neutral ino_l 


tau"- 


tau"+ 


0. 


. 00000422 


0. 


.02157388 


Total width : 






0. 


.01956317 







Neutralino_3 
Selectron~-_l e~+ 
Selectron~+_l e~- 
Selectron~-_2 e~+ 
Selectron~+_2 e~- 
Smuon~-_l mu~+ 
Smuon~+_l mu~- 
Smuon~-_2 mu"+ 
Smuon~+_2 mu~- 
Stau~- 1 tau~+ 



0.00246761 
0.00246761 
0.00110911 
0.00110911 
0.00249046 
0.00249046 
0.00115204 
0.00115204 
0.01021462 



0.12381862 
0.12381862 
0.05565227 
0.05565227 
0.12496504 
0.12496504 
0.05780683 
0.05780683 
0.51254502 
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Stau~+_1 tau~- 


0. 


.01021462 


0. 


.51254502 


Stau~-_2 tau~+ 


0. 


,01218195 


0. 


.61126086 


Stau~+_2 tau~- 


0. 


.01218195 


0. 


.61126086 


e-sneutrino nu_e~* 


0. 


.00612932 


0. 


.30755439 


e-sneutrino~* nu_e 


0. 


.00612932 


0. 


.30755439 


mu-sneutrino nu_mu~* 


0. 


.00612947 


0. 


.30756211 


mu-sneutrino ~* nu_mu 


0, 


.00612947 


0. 


.30756211 


tau-sneutrino nu_tau~* 





.00617278 


0. 


.30973532 


tau-sneutrino ~* nu_tau 


0. 


.00617278 


0. 


.30973532 


Chargino~+_l W~- 


0. 


. 59455684 


29. 


.83342506 


Chargino~-_l W~+ 


0. 


. 59455684 


29. 


.83342506 


Neutralino_l Z 


0. 


.22190910 


11. 


. 13486222 


Neutralino_2 Z 


0. 


.41946063 


21. 


. 04752045 


Neutral ino_l hO 


0. 


.04208191 


2. 


.11156877 


Neutralino_2 hO 


0. 


.02422138 


1. 


.21537033 


Neutralino_2 photon 


0. 


.00002024 


0. 


.00101568 


neutralino_l b b~* 


0. 


.00000772 


0. 


.00038719 


neutralino_2 b b~* 


0, 


.00000522 


0, 


.00026192 


Total width : 


1 


.99292182 






Neutral ino_4 










Selectron~-_l e~+ 


0. 


.01008587 


0. 


.38204392 


Selectron~+_l e~- 


0. 


.01008587 


0. 


.38204392 


Selectron~-_2 e~+ 


0. 


.02480502 


0. 


.93959259 


Selectron~+_2 e~- 


0. 


.02480502 


0. 


.93959259 


Smuon~-_l mu~+ 


0. 


.01007017 


0. 


.38144944 


Smuon~+_l mu~- 


0, 


.01007017 


0. 


.38144944 


Smuon~-_2 mu~+ 


0, 


.02486868 


0. 


.94200385 


Smuon~+_2 mu"- 


0, 


.02486868 


0. 


.94200385 


Stau"-_1 tau"+ 


0. 


.00789319 


0. 


.29898723 


Stau~+_1 tau~- 


0. 


.00789319 


0. 


.29898723 


Stau"-_2 tau"+ 


0. 


.04123849 


1. 


.56207842 


Stau"+_2 tau~- 


0. 


.04123849 


1. 


.56207842 


e-sneutrino nu_e~* 


0. 


.06494620 


2. 


.46010577 


e-sneutrino~* nu_e 


0. 


.06494620 


2. 


.46010577 


mu-sneutrino nu_mu~* 


0. 


. 06494764 


2. 


.46016018 


mu-sneutrino ~* nu_mu 


0. 


. 06494764 


2. 


.46016018 


tau-sneutrino nu_tau~* 


0, 


.06535235 


2. 


.47549054 


tau-sneutrino"* nu_tau 





.06535235 


2. 


.47549054 


Chargino"+_l W"- 


0. 


.67613622 


25, 


.61145377 


Chargino~-_l W"+ 





.67613622 


25. 


.61145377 


Neutralino_l Z 


0. 


.05545958 


2, 


.10076061 


Neutralino_2 Z 


0. 


. 04990974 


1. 


.89053781 


Neutral ino_l hO 


0. 


. 18282528 


6. 


.92526289 


Neutralino_2 hO 


0. 


.37096704 


14. 


.05190972 


Neutralino_2 photon 


0. 


.00002906 


0. 


.00110069 
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chargino~+_l d u~* 
chargino~-_l d~* u 
chargino~+_l s c~* 
chargino~-_l s~* c 
chargino~+_l b t~* 
chargino~-_l b~* t 
neutralino_l b b~* 
neutralino_2 u u~* 
neutralino_2 c c~* 
neutralino_2 d d~* 
neutralino_2 s s~* 
neutralino_2 b b~* 
Total width : 



0.00000417 
0.00000417 
0.00000418 
0.00000418 
0.00001489 
0.00001489 
0.00000829 
0.00000460 
0.00000459 
0.00000530 
0.00000530 
0.00001293 
2.63997596 



0.00015812 
0.00015812 
0.00015820 
0.00015820 
0.00056401 
0.00056401 
0.00031410 
0.00017408 
0.00017392 
0.00020077 
0.00020085 
0.00048988 



Gluino 
Sup_l u~* 
Sup_l~* u 
Sup_2 u~* 
Sup_2~* u 
S-charm_l c~* 
S-charm_l~* c 
S-charm_2 c~* 
S-charm_2~* c 
Stop.l t~* 
Stop.l"* t 
Sdown.l d~* 
Sdown_l~* d 
Sdown_2 d~* 
Sdown_2~* d 
S-strange_l s' 
S-strange_l~* 
S-strange_2 
S-strange_2~* s 
Sbottom.l b~* 
Sbottom_l~* b 
Sbottom_2 b~* 
Sbottom_2~* b 
Stop.l c~* 
Stop.l"* c 
neutral ino_l gluon 
neutralino_2 gluon 
neutralino_3 gluon 
neutral ino_4 gluon 
neutral ino_l t t~* 
neutralino_2 t t~* 
chargino"+_l t"* b 



0.22870900 
0.22870900 
0.11126815 
0.11126815 
0.22850951 
0.22850951 
0. 11134187 
0. 11134187 
0.23793220 
0.23793220 
0.23074812 
0.23074812 
0.08317483 
0.08317483 
0.23077130 
0.23077130 
0.08317883 
0.08317883 
0.51851488 
0.51851488 
0.24263650 
0.24263650 
0.00530170 
0.00530170 
0.00002680 
0.00022135 
0.00033270 
0.00038117 
0.00008962 
0.00009374 
0.00066796 



4.94258582 
4.94258582 
2.40459439 
2.40459439 
4.93827473 
4.93827473 
2.40618760 
2.40618760 
5. 14190657 
5.14190657 
4.98665276 
4.98665276 
1.79747505 
1.79747505 
4.98715370 
4.98715370 
1.79756163 
1.79756163 
11.20552447 
11.20552447 
5 . 24357028 
5 . 24357028 
0. 11457405 
0. 11457405 
0.00057923 
0.00478362 
0.00718998 
0.00823738 
0.00193684 
0.00202577 
0.01443507 
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chargino~-_l t b~* 
chargino~+_2 t~* b 
chargino"-_2 t b~* 
Total width : 



0.00066796 
0.00032981 
0.00032981 
4.62731469 



0.01443507 
0.00712747 
0.00712747 



hO 

muons 
taus 
d-quark 
s-quark 
b-quark 
c-quark 
Total width 



0.00000101 
0.00028643 
0.00000000 
0.00000104 
0.00180832 
0.00009191 
0.00218871 



0.04633093 
13.08657852 

0.00016098 

. 04730787 
82.62040292 

4. 19919143 



HO 



muons 




0, 


.00027932 


0. 


,03671913 


taus 




0. 


.07900846 


10. 


.38638120 


d-quark 




0, 


.00000098 


0. 


00012915 


s-quark 




0, 


.00028872 


0. 


03795547 


b-quark 







.50606084 


66. 


,52630492 


c-quark 




0, 


.00000440 


0. 


00057844 


t-quark 




0. 


.03284785 


4. 


,31814858 


Selectron 1 1 




0. 


.00038588 


0. 


,05072782 


Smuon 1 1 




0, 


.00039320 


0. 


,05168917 


Smuon 1 2 







.00001885 


0. 


, 00247807 


Smuon 2 1 




0, 


.00001885 


0. 


, 00247807 


Stau 1 1 




0, 


.00420030 


0. 


,55216824 


Stau 1 2 




0, 


.00406411 


0. 


,53426431 


Stau 2 1 







,00406411 


0. 


.53426431 


e-Sneutrino 




0, 


.00071373 


0. 


.09382650 


mu-Sneutrino 




0, 


.00071386 


0. 


. 09384330 


tau-Sneutrino 




0. 


.00074895 


0. 


.09845649 


neutralino_l neutralino. 


.1 


0. 


.01640703 


2. 


. 15685379 


neutralino_l neutralino. 


.2 


0. 


.04741210 


6. 


.23275196 


neutralino_2 neutralino. 


.2 


0. 


.01321269 


1. 


.73692880 


chargino~+_l chargino"-. 


.1 


0. 


.03572858 


4. 


.69684775 


Z Z 




0, 


,00145341 


0. 


.19106428 


w+ w- 




0. 


.00311215 


0. 


.40912109 


hO hO 




0. 


.00955444 


1. 


,25601818 


Total width : 




0. 


.76069285 







AO 

muons 0.00027991 0.02394510 

taus 0.07918234 6.77364982 
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s-quark 








0, 


.00028933 


0. 


.02475108 


b-quark 








0. 


.50718753 


43, 


.38733581 


c-quark 








0, 


.00000337 


0. 


.00028798 


t-quark 








0, 


. 10446570 


8, 


.93651435 


Smuon 1 2 








0. 


.00001955 


0, 


.00167248 


Smuon 2 1 








0. 


.00001955 


0. 


.00167248 


Stau 1 2 








0, 


.00569243 





.48695874 


Stau 2 1 








0, 


.00569243 


0. 


.48695874 


neutral ino. 


.1 


neutralino. 


.1 


0. 


.02519311 


2 


.15514326 


neutral ino. 


.1 


neutralino. 


.2 


0. 


.10809913 


9 


.24733561 


neutral ino. 


_2 


neutralino. 


.2 


0. 


.09248237 


7, 


.91140044 


chargino~+_ 


.1 


chargino"-. 


.1 


0. 


.23769243 


20 


.33338946 


hO Z 








0, 


.00267578 


0. 


.22889978 


Total width 








1. 


. 16897596 







H~+ 



muon neutrino 





.00028573 


0. 


,04266509 


tau neutrino 


0. 


.08082739 


12. 


,06922654 


d-quark u-quark 


0. 


.00000088 


0. 


.00013177 


s-quark c-quark 


0. 


.00026233 


0. 


.03917183 


b-quark t-quark 


0. 


.42729169 


63. 


.80362064 


Selectron_2 Sneutrino 


0. 


.00051866 


0. 


. 07744696 


Smuon_l Sneutrino 


0. 


.00005543 


0. 


.00827732 


Smuon_2 Sneutrino 


0. 


.00051472 


0. 


.07685861 


Stau_l Sneutrino 


0. 


.01489058 


2. 


.22347616 


Stau_2 Sneutrino 


0, 


.00000817 


0, 


.00121960 


chargino_l neutralino_l 


0, 


. 14117172 


21. 


.07990203 


chargino_l neutralino_2 


0, 


.00083132 


0. 


.12413279 


hO W 





.00303956 


0. 


.45386947 


Total width : 


0. 


.66969819 







Total cross sections in fb for: 
E_cms : 500.0000000 GeV 

Degree of polarization: P_e- = 0.000000 P_e+ = 0.000000 
Intial state radiation is included 

u-squarks : kinematically not possible 

c-squarks : kinematically not possible 

t-squarks : kinematically not possible 



d-squarks : kinematically not possible 
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s-squarks : kinematically not possible 
b-squarks : kinematically not possible 



selectrons 



1 1 


283. 


.6216345 


fb 


1 2 


78. 


.0199000 


fb 


2 2 


45. 


.7859068 


fb 


smuons 








1 1 


54. 


.7929652 


fb 


1 2 


0. 


,0059693 


fb 


2 2 


19. 


,0224584 


fb 


staus 








1 1 


59. 


.6056340 


fb 


1 2 


1. 


.2593673 


fb 


2 2 


15. 


.9954658 


fb 


e-sneutrino 








453. 


.2694763 


fb 


mu-sneutrino 








13, 


.6239191 


fb 


tau-sneutrino 








13. 


.9240765 


fb 


neutralinos 






1 1 


275. 


.7338081 


fb 


1 2 


65, 


, 7329474 


fb 


1 3 


7. 


.2644701 


fb 


1 4 


0. 


.9005054 


fb 


2 2 


69. 


.0994560 


fb 


charginos 








1 1 


161. 


,0778560 


fb 



h~0 Z 



24.7411048 fb 
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H~0 Z : cross section below 1 . 0000000E-03 fb 
h~0 A~0 : kinematically not possible 
H~0 A~0 : kinematically not possible 
H~+ H~- : kinematically not possible 
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